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ABSTRACT 

A  finite  element  code  is  developed  for  analysis  and 
design  of  three  dimensional  truss  and  frame  structures. 
Structures  are  designed  for  minimum  weight  subject  to 
constraints  on:   member  stresses,  Euler  buckling,  shell 
buckling,  joint  displacements  and  system  natural  fre- 
quencies .   Structures  are  optimized  with  respect  to  member 
size  and  structure  configuration. 

The  finite  element  code  may  be  used  for  analysis  only, 
or  may  be  coupled  to  an  optimizer  of  the  user's  choice. 
The  displacement  method  is  used  for  static  analysis,  and 
structure  natural  frequencies  are  calculated  via  the  sub- 
space  iteration  method. 

Flexibility  is  provided  for  expansion  to  other  than 
tubular  frame  elements,  and  provisions  are  made  for  the 
future  growth  to  panel  and  other  types  of  structural 
elements . 

Documentation  is  provided  to  facilitate  use  of  the 
code.   A  User's  manual  is  presented  with  examples  and  re- 
sults .   An  explanation  of  how  this  code  may  be  coupled  to 
an  optimizer  is  also  provided. 
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I.   INTRODUCTION 

The  task  of  the  Engineer  is  to  provide  the  best  solution 
to  the  customer's  problem.   Usually  the  "best"  solution  is 
the  one  that  does  the  job  with  an  adequate  margin  of  safety, 
is  aesthetically  pleasing,  and  is  economically  feasible. 
The  solution  may  be  reached  through  various  means ,  but  effi- 
cient use  of  design  tools,  as  well  as  efficient  use  of 
materials  is  important  since  both  add  to  the  overall  cost  of 
the  product.   Design  optimization  is  one  method  that  can  be 
used  to  maximize  the  efficiency  of  a  structure  by  minimizing 
its  weight  and,  presumably,  cost. 

Optimization  of  structures  has  had  continuing  changes 
since  its  development  in  the  early  1960  '  s  with  an  active 
area  of  research  being  elastic  structures.   The  main  goal  is 
to  design  structural  systems  that  efficiently  perform  speci- 
fied purposes.   Since  most  physical  structures  can  be 
modeled  by  the  finite  element  method,  a  computer  program 
can  be  written  to  perform  the  necessary  calculations  to 
solve  the  problem. 

The  ourpose  of  this  research  was  to  develop  a  finite 
element  code  that  could  be  used  to  analyze  a  combined  truss/ 
frame  structure  and  could  be  easily  coupled  to  an  optimizer; 
thus  croviding  a  useful  tool  for  designing  such  structures 
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as,  for  example,  ships'  masts.   This  code  expands  the  pre- 
vious work  by  Fitzgerald  [Ref .  1]  on  truss  structures  to 
the  more  general  six  degree-of-freedom  truss/frame  case. 

The  design  problem  considered  in  this  study  is  the  opti- 
mization of  three-dimensional  statically  indeterminant  com- 
bined truss-frame  structures  under  varying  load  conditions 
and  subject  to  a  variety  of  constraints.   The  objective  is 
to  minimize  the  weight  of  the  structure  where  the  design 
variables  are  member  sizes  and  joint  coordinates.   Constraints 
include:   maximum  normal  stress;  maximum  joint  displacements; 
minimum  structure  natural  frequencies;  Euler  buckling;  and 
in  rhe  case  of  the  tubular  frame  elements ,  local  or  shell 
buckling . 

In  the  present  code,  all  gradient  information  is  calcu- 
lated by  the  finite  difference  method.   Modification  of  the 
code  to  permit  calculation  of  gradients  analytically  has 
been  identified  as  a  necessary  future  extension. 

This  document  describes  the  use  and  capabilities  of  the 
finite  element  computer  code  to  be  coupled  to  an  optimizer. 
The  user's  manual  presented  in  Chapter  V  contains  a  simple 
design  example  in  which  the  program  is  coupled  to  the  CONMIN 
optimization  code  [Ref.  23.   Additionally,  guidelines  for 
coupling  the  code  to  an  optimizer  of  the  user's  choice  are 
presented . 
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Several  examples  demonstrating  the  program  under  a 
variety  of  conditions  are  presented.   Conclusions  and  recom- 
mendations for  future  work  are  given. 
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II.       ANALYSIS 

A.  INTRODUCTION 

When  the  finite  element  method  of  analysis  is  used  to 
design  optimization,  two  objectives  must  be  kept  in  mind. 
First,  the  number  of  analyses  for  the  structure  should  be 
kept  to  a  minimum.   Second,  the  amount  of  gradient  informa- 
tion required  during  the  design  process  should  be  minimized 
to  shorten  run  times  and  minimize  computer  storage 
requirements . 

B.  STATIC  ANALYSIS 

Initial  formulation  of  the  problem  must  include  approxi- 
mate member  areas  in  the  case  of  truss  elements,  and  for 
frame  elements,  characteristic  dimensions  (for  tubular  members 
mean  diameter  and   wall  thickness);  material  properties  (which 
may  be  different  for  each  member) ;  a  set  or  sets  of  external 
loads;  any  non-structural  attached  masses;  and  specified 
joint  support  conditions . 

The  analysis  for  the  stresses  and  deflections  at  the 
joints  must  satisfy  the  conditions  of  equilibrium  of  forces 
at  the  nodes  and  geometric  conditions  of  compatibility  of 
deformation.   In  this  analysis  the  structure  is  assumed  to 
behave  in  a  linearly  elastic  fashion.   The  weight  of  an 
individual  member  is  not  inherently  included  as  part  of  the 
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specified  load  conditions;  but,  as  an  option,  half  of  the 
weight  of  each  member  may  be  applied  at  the  member's  end- 
points  as  an  additional  load  in  the  negative  Y-direction. 

For  this  analysis,  the  following  assumptions  are  made: 
truss  and  frame  members  are  treated  as  discrete  entities; 
truss  elements  have  three  translational  degrees  of  freedom 
at  each  node  and  are  treated  as  pin-connected;  frame  elements 
have  three  translational  and  three  rotational  degrees  of 
freedom  at  each  node  and  are  treated  as  fixed-fixed  beams ; 
and  loads  and  reactions  are  applied  at  the  joints  as  shown 
in  Figure  2.1. 

The  Displacement  (Stiffness)  method  for  finite  element 
analysis  [Ref .  3]  is  utilized  where 

Ku  =  F  (3.1) 

and  where  K  is  the  global  stiffness  matrix,  J  is  the  vector 
or  vectors  of  aoplied  loads,  and  u  is  the  vector  or  vectors 
of  displacements.   The  method  used  herein  is  an  extension  of 
that  described  by  Felix  and  Vanderplaats  [Ref.  4].   3y 
applying  the  constitutive  stress-displacement  relationships , 
stresses  in  the  elements  may  be  recovered. 

C.   DYNAMIC  ANALYSIS 

When  constraints  on  the  system's  natural  frequencies 
are  to  be  considered,  the  design  process  requires  the  solu- 
tion of  an  eigenvalue  problem.   The  sub-space  iteration 
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Figure  2.1   FORCE  ORIENTATION  CONVENTION  FOR  AN  ELEMENT 
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method  of  Bathe  and  Wilson  [Ref .  5]  is  used  to  solve  for  the 
desired  number  of  lowest  eigenvalues  and  the  associated 
eigenvectors .   This  method  is  reasonably  efficient  for  a 
small  number  of  lowest  frequencies  for  a  large  problem,  and 
is  well  suited  for  re-analysis  when  small  changes  are  made 
in  the  design. 

D.   GRADIENTS 

Gradients  are  currently  calculated  with  respect  to  member 
sizes  and  coordinates  only  by  finite  difference  techniques. 
Inclusion  of  the  capability  to  calculate  gradients  analytic- 
ally is  identified  as  a  needed  extension  to  this  work. 
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III.   OPTIMIZATION 

A.   INTRODUCTION 

The  primary  objective  of  structural  optimization  is  to 
design  systems  that  efficiently  perform  specified  purposes. 
Selection  of  a  specific  optimizing  algorithm  must  take  into 
account  the  following:   1)  the  structure  should  be  analyzed 
as  few  times  as  possible;  2)  the  algorithm  should  minimize 
the  amount  of  gradient  information  required;  3)  the  algorithm 
should  provide  reasonable  assurance  that  an  optimum  or  near- 
optimum  design  will  be  reached. 

3.   GENERAL  FORMULATION 

The  general  statement  of  inequality  constrained  minimiza- 
tion is  as  follows: 
Minimize 

F(X)  («*.l) 

Subject  to: 

G-(X)   £   0  3=1, m  (4.2) 

xl    <    x.  <  xu  i=l,n  (4.3) 

i  —   i  —   i 

where  F(X)  is  the  objective  function  to  be  minimized.   The 

functions  G.(X)  are  the  set  of  inequality  constraints  to 

j  ~ 

be  met.   The  vector  X  contains  the  design  variables. 
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The  inequality  constraints,  G-(X)  <_  0   j=l,m,  must  be 

satisfied  for  the  design  to  be  accepted  as  feasible.   Sice 

I  u 

constraints  ,  x^  and  x.  ,  are  lower  and  upper  bounds  on  the 

design  variables.   The  objective  function  must  be  minimized 

as  much  as  possible  while  still  satisfying  all  inequality 

constraints.   In  the  case  that  it  is  not  possible  to  satisfy 

all  constraints,  the  most  nearly  feasible  solution  must  be 

found.   Felix  and  Vanderplaats  [Ref.  4]  is  an  excellent 

source  for  the  basic  structural  design  formulation. 

C.   DESIGN  VARIABLES 

The  vector  X  contains  the  design  variables;  in  this  case, 
member  characteristic  dimensions  for  frame  elements,  cross- 
sectional  areas  for  truss  elements ,  and  spacial  coordinates 
of  the  joints.   The  user  may  elect  to  optimize  the  structure 
weight  with  respect  to  any  combination  of  the  above  variables 

For  truss  elements,  where  normal  stress  is  dependent  only 
upon  the  magnitude  of  the  cross-sectional  area  of  an  element 
and  not  the  distribution  of  that  area,  it  is  sufficient  to 
use  the  area  as  the  design  variable  as  long  as  the  Euler 
buckling  stress  can  be  related  to  the  cross-sectional  area. 
For  frame  elements  where  stresses  are  dependent  on  area 
distribution,  more  flexibility  is  allowed  by  varying  two 
characteristic  dimensions  and  allowing  the  code  to  calculate 
the  section  oroDerties  from  these . 
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In  the  case  of  tubular  elements,  the  characteristic 
dimensions  are  mean  diameter  and  wall  thickness;  from  which 
area,  bending  and  polar  moments,  and  maximum  radial  dimen- 
sions are  generated. 

The  optimum  geometry  problem  requires  that  the  joint 
coordinates  be  design  variables.   The  x,  y,  and  z  coordinates 
of  a  joint  are  treated  as  separate  design  variables. 

In  many  cases  it  is  desirable  to  link  a  set  of  design 
variable  together  to  preserve  symmetry,  limit  the  number  of 
variables  to  be  solved,  or  to  limit  the  number  of  unique 
elements  to  be  manufactured.   The  code  has  provisions  for 
design  variable  linking  by  which  two  or  more  variables  may 
be  linked  in  equality  or  some  fixed  ratio. 

D.   OBJECTIVE  FUNCTION 

The  objective  function  under  consideration  is  weight 


NE 
F(X)   =     S    P.A-L.  (4.4) 

i  =  l   " 


where  p.  is  the  weight  density  (in  consistent  units)  of  the 

J- 

material  of  the  ith  element,  A.  is  the  cross-section  area 
of  the  ith  element,  and  L-  is  the  length  of  the  ith  element 
and  NE  is  the  number  of  elements  in  the  structure. 

E .   CONSTRAINTS 

This  code  is  designed  to  accommodate  constraints  on  maxi 
mum  normal  stress,  Euler  buckling,  local  or  shell  buckling, 
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maximum  joint  displacements,  and  minimum  natural  fre:.^:.:.-: 
of  the  structure.   All  constraint  values  are  normalized. 
1 .   Stress : 

For  truss  elements,  maximum  normal  stress   J,,   is 
calculated  directly  from  the  element  tensile  forces  .   7or 
frame  elements;  tensile  stress  o*     ,  maximum  bending  end 
shear  stresses   J«  and  o„   are  calculated.   Maximum  and 
minimum  normal  stresses  are 


.  +  /(a.    -    '^a^ 


o  -    max  magnitude  ::  J,  +  /(a,    *  4a^)       (4.5) 

max         H 


=  min  magnitude  of  a   +  /(c    +  imj   )      (4.5) 


mm      = A  - 


The  upper  and  lower  bounds  on  stress  may  be  different  for 
each  member,  but  are  taken  to  be  the  same  for  every  loading 
condition . 

2 .   Local  cr  Thin  Shell  Suckling 

The  stress  at  which  local  or  thin  shell  buckling 
occurs  is  given  by: 

a.      0.4E-  (4.7) 


m-  /  t« 


where    the    subscript    i    ccrresoonds    to    the    tether    nutter, 
is    Young's    modulus.    1„      is    the    element's    teen    diameter , 

t-    is    the    element's    wall    thickness. 
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F.   DESIGN  VARIABLE  BOUNDS 

Side  constraints  are  imposed  on  the  design  variables 
as  : 

CD*    £   CD.    £   CD^  (4.3) 

I  u 

CD  is  the  characteristic  dimension  and  CD-  and  CDT  are 

the  minimum  and  maximum  allowable  characteristic  dimensions 
of  the  ith  element,  and  are  taken  to  be  the  same  for  all 
load  cases.   If,  as  is  the  case  in  a  tubular  element,  geo- 
metry dictates  that  some  relationship  between  the  design 
variables  cannot  be  exceeded,  i.e.,  thickness  cannot  exceed 
the  mean  diameter,  the  user  must  arrange  the  bounds  to  pre- 
clude such  an  occurrence . 


OPTIMUM  GEOMETRY 


Joint  coordinates  are  treated  as  design  variables  with 
a  separate  design  variable  for  each  coordinate  direction 
of  a  node.   Coordinate  design  variables  may  also  be  linked 
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IV.   PROGRAM  FEATURES 

A.   INTRODUCTION 

Computer  codes  each  have  their  own  features  and  formats 
with  which  the  user  must  become  familiar  if  he  is  to  use 
the  code  easily  and  efficiently.   The  SADX  code  developed 
in  the  course  of  this  research  has  been  designed  with 
various  options  which  are  explained  in  this  chapter. 
Chapter  V  contains  a  User's  Manual  with  sample  data  for  a 
typical  problem  that  might  be  solved  with  this  code:   a 
truss-braced  cantilever  beam.   This  problem  along  with 
other  numerical  examples  will  be  presented  in  detail  with 
results  in  Chapter  VI. 

The  SADX  code  was  written  to  be  used  as  a  stand  alone 
analysis  program  or  as  an  analysis  code  that  could  be  easily 
be  coupled  to  an  optimizer  (of  the  user's  choice)  through 
simple  modifications  to  the  main  driver  program. 

3 .   FEATURES 

When  the  user  supplies  member  areas,  section  types, 
characteristic  dimensions,  material  data,  connectivities 
and  joint  coordinates,  along  with  a  set  of  program  control 
parameters;  the  analysis  mode  will  calculate  the  weight  of 
the  structure.   The  addition  of  one  or  more  sets  of  loading 
parameters  will  result  in  the  calculation  of  resultant  joint 
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displacements,  member  stresses,  and/or  forces  for  each  load 
case  along  with  the  desired  number  of  structure  natural 
frequencies  and  modes.   Design  variables  may  be  chosen  as 
truss  element  areas,  frame  element  characteristic  dimen- 
sions, and  joint  coordinates.   When  coupled  to  an  optimizer, 
the  code  will  minimize  the  weight  of  the  structure  and  print 
the  final  optimization  information.   The  code  is  designed  to 
be  as  simple  to  use  as  possible  while  retaining  the  flexi- 
bility for  use  on  a  variety  of  problems . 

The  code's  modular  construction  allows  the  user  to  use 
frame  element  cross-sections  other  than  tubular  elements. 
This  is  done  by  reading  in  two  characteristic  dimensions  for 
each  frame  element  along  with  a  section  type  identifier. 
Subroutine  SADX35  is  called  to  calculate  the  section  proper- 
ties, area,  maximum  radial  dimensions,  and  bending  and  polar 
moments  of  inertia.   The  user  may  augment  this  subroutine  to 
calculate  section  properties  for  whatever  section  type  he 
may  choose  to  work  with. 

Many  computer  codes  require  the  definition  of  an  auxil- 
iary node  to  orient  the  principle  axis  of  a  non-axially 
symmetric  element.   In  this  code  it  is  assumed  that  the 
element's  local  z-axis  is  parallel  to  the  global  x-z  plane. 

Pseudo-dynamic  storage  is  used  to  allow  storage  of  most 
data  in  one  integer  and  one  real  array  for  more  efficient 
use  of  storage . 
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Various  print  control  options  are  available  to  tailor 
the  printed  output  to  match  the  user's  desires. 

Design  variable  linking  is  available  to  allow  elements 
and  joints  to  be  grouped  to  maintain  symmetry,  limit  the 
number  of  independent  design  variables,  or  reduce  the 
variety  of  member  sizes  generated. 

Optimization  may  be  performed  with  respect  to  member 
size,  with  respect  to  structure  geometry,  or  with  respect  to 
both. 

Structures  may  be  comprised  of  truss  elements,  frame 
elements,  or  a  mix  of  the  two  types. 

Structures  may  be  optimized  for  multiple  load  cases  with 
constraints  imposed  upon  any  combination  of  maximum  normal 
stress,  maximum  joint  displacements  and  rotations,  Euler 
and  local  buckling,  and  a  specified  number  of  minimum  natural 
frequencies  in  free  vibration.   Separate  displacement  con- 
straints may  be  imposed  for  each  load  case  .   3oth  consis- 
tent and  lumped  mass  options  are  coded.   Either  forces  or 
stresses  or  both  can  be  output.   The  user  may  decide  whether 
to  include  the  structure's  weight  and  the  fixed  masses  as 
loads  applied  to  the  structure. 

C .   EXAMPLE 

The  following  example  of  the  truss-braced  cantilever 
beam  presented  in  Tables  (I-V)  demonstrates  some  of  the 
options  available  in  the  code. 
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TABLE    I 
TRUSS -BRACED   CANTILEVER   BEAM:       INPUT    CONTROL    PARAMETERS 


INPUT    PARAMETERS    FOR    STRUCTURAL    ANALYSIS  AND    DESIGN 
ROUTINE,    "SADX" 

TRUSS-ERACED    CANTILEVER     BEAM:    EXAMPLE    1 

CONTROL    PARAMETERS 

TOTAL    NUMBER    OF    ELEMENTS,  NE=  4 

TOTAL    NUMBER    OF    BAR    ELEMENTS,  NEB  =  2 

TOTAL    NUMBER    OF    FRAME    ELEMENTS,  NEF=  2 

TOTAL    NUMBER    OF    JOINTS,  NJ=  5 

JOINT    MAXIMUM    DEGREES    OF     FREEDOM,     JCM=  6 

JOINT    CONSTRAINT    VARIABLE,  NCJ=  3 

NUMBER    OF    MATERIAL    TYPES,  NMT=  2 

NUMBER    OF    LOAD    CONDITIONS,  NLC=  2 

NUMBER    OF    EIGENVALUES,  NEIG=  1 

NO.    OF    EIGENVALUES    CALCULATED,       NEIG1=  2 

NUMBER    OF    FIXED    MASSES,  NFSASS=  1 

EULE?    BUCKLING    CONSTRAINT    ID,         NE03C=  1 

LOCAL    BUCKLING    CONSTRAINT    ID,         L3UCK=  1 

NO.    OF    DI5PL.    CONSTRAINTS,  ND5?LC=  2 

NO.    OF    FREQ.     CONSTRAINTS,  NFR2Q=  1 

LUMPED    MASS    OPTION,  L1ASS=  0 

FORCE/STRESS    PRINT    OPTION,  NSIRES=  2 

OPTIMUM    SIZS/BOTH/GEOM    OPTION       Id7CLC=  2 

STRUCTURE    WEIGHT    AS    LOADS    OPTION    NSTW=  1 

FIXED    MASSES    AS    LOADS    OPTION  NFM'w  =  1 

ACCELERATION    DUE    TO   GRAVITY,     GRA7  =    0.38640E+03 

EIGENVALUE    CONVERGENCE    TOLERANCE,     EP3EIG  =    0. 100002-03 
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TABLE    II 
EXAMPLE    1:       JOINT    COORDINATE   AND   MEMBER   INPUT   DATA 


JCINT 

CCCR0IN*TES 

JCINT 

X 

Y 

I 

1 

o.c 

0.0 

0 

.0 

2 

O.1C00OE+03 

0.0 

0 

•  0 

3 

0.2CCCOE+C3 

0.0 

0 

.0 

4 

0.0 

C.  150O0E+0  2 

0 

.0 

5 

0.0 

0.0 

0 

.500  00E+02 

CCORDINATE 

variable: 

CESIGN 

VARIABLE 

MULTIPLIER 

JCINT 

X      Y       2 

X 

Y 

z 

1 

C      0       C 

0.0 

0.0 

0.0 

2 

C      0      C 

0.0 

0.0 

0.0 

3 

C     0      c 

0.0 

0.0 

0.0 

4 

C     1      c 

0.1CC0E+01       0.1000E+ 

01 

0.1000E+01 

5 

C      0       2 

C.10COE+01       0.100CE+ 

01 

O.IQOOE+Oi 

ELEMENT    INFOFPATICN    PCP    EAR    ELEMENTS 

ELEMENT-JOINT    RELATIONSHIPS 

LNG    N0DE1    N0DE2    MATL    0VAR1  AREA 

12  4  11         0.  3000E+C1 

2  2  5  12         0.3000E+C1 


LENGTH 
0.  1303E+03 
0.  1118E+03 


ELEMENT    INFORPATICN    FCP    FRAME    ELEMENTS 

ELEMENT-JOINT    RELATICNSHI  PS         ELEMENT    PROPERTIES 

ELEMENT    NUMBEF  3 

NCDE1       NGDE2       MATL      1ST    OESVAR    2ND    DESVAR 
12  2  3  4 

AREA  LENGTH  Z-MOMENT         Y-MOMENT 

0.157E+C2       0.100E+02      0.511E+02      0.511E+02 


ELEMENT    NUMBEF  A 

NCDE1       N0DE2       PATL      1ST    OESVAR    2ND    CESVAR 
2  3  2  3  5 

AREA  LENGTH  Z-MOMENT         Y-MOMENT 

0.157E-K2      0.1CCE+C2      C.511E+02      0.511E+02 


ZMAX 
0.300E+01 


ZMAX 
0.3005+01 


YMAX 

0.300EH 


YMAX 

0.300E+1 
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TABLE  III 
EXAMPLE  1:   SAMPLE  DISPLACEMENT  AND  FORCE/STRESS  OUTPUT 


»c.  e»  (j»  eliments  • 


INC 

I 

1 

10*0    COO 
t 
2 

TENSILE    FORCE     IN    X 

-C.73C7F.O* 

0.  1S*9E»0« 

N»»    STRESS 
-O.tf28F«03 

LNC 
1 

1.0*0   CCKO 
2 

TENSILE    FORCE    IN    X 

0.7*20E»0* 

-0.2<.136*OS 

M»X    STRESS 

O.87*o£»03 
-0).*016E»0* 

NC.    OF 

EH*Nf    ELEMENTS    • 

2 

IllNERT 

LO*C    CONO 

NCOE 

TENSKE    FCRCE       SHEAR    FORCE 
IN    X    OIREX               IN    r    0IR6X 

SMEAR    FORCE 
IN    2    QIREX 

TCRSICN    MOMENT 
AdOUT    X    AXIS 

8EN0ING   NCHENT 
A8CUT    Y    tXIS 

RENO  I  *G    MCMENT 
A8CUT    I    AXIS 

MAX 

STRESS 

I 

3 

1 
1 

ICk 
HIGH 

-0.1T57E«0* 
0.1757E»04 

0.10666*0* 
-0. 10666*06 

-0.33956*01 
0.33956*01 

0.0 

0.0 

0.3395E*C3 
-0.97666-03 

0.35156.05         -c. 
0.  71*56*05           0. 

«6?:e*o* 

19C3E*05 

3 
3 

2 
2 

ICk 
H1CH 

0. 66*76 »0* 
-O.66*7E»0* 

-0.»18TE*03 

3.-J37E.03 

-0.1*826*0* 

o.i*826«o* 

0.0 

0.0 

0.*816E*C5 
0.10C06*C6 

-C.13i?E»05           0. 
■0.23556*05        -0. 

1**76*05 
27«GE*0S 

element 

LO»C   CCNO 

kCBC 

TENSILE   FCRCE      SHEAR    FORCE 
IN   X   OIREX             IN  1    OIREX 

SMEAR    FORCE 
IN    2    SIREX 

TOHStCN    FOMENT 
A8CUT    X    AXIS 

5EK0ING    MOMENT 
ABCUT    Y    AXIS 

8SN0IK    MCM6NT 
ABOUT    ;    AXIS 

MAX 

STRESS 

» 
* 

1 

1 

ICk 

HIGH 

-0.*88?E-03 
0.*882E-O3 

-0.71*56*03 

0.71*5E«O3 

-0.*768E-05 
0.4768 E-OS 

0.0 
0.0 

0.13  73E-03 
0.732*6-03 

•0.71*56*05         -0 
0.6250E-01           0. 

2*726*05 
2271E-01 

*  2                    LC»               0. 39066-02 

*  I                    HIGH         -0.39066-02 
EIGENVALUES    1IC    EIGENVECTORS 

0.26556*03 

-O.J  S3  56*03 

C.10OOE*0* 

-0. .000  6*0* 

0.0 

o.o 

-0.10006*06 

-O.W75E*GO 

0.2S55E-O5            0 
-0. 625C6-01         -0. 

.36006*05 
6915E-01 

HCOE    MF8ER            1 
EFECUERCT      •    C.*ie8SE»02 
C.2631'iE»03 
EIGENVALUE    •    C.<927CE»C5 
EICENVICTOR       C6GREE    CF    F 
JCIST                 1 

1  0.0                               0. 

2  3.0                               C. 

3  -C.2?«:5E-fll       -C. 
*         -0.50717E-01       -C. 
5          -C. 23*356-01       -0. 

CM 

SE-S 

>EEMN 

2 
1                            0. 
}                            0. 

37SME-01       -0. 
IC7456-01          0. 
171S5E-01          0. 

3 

0 

J 

74S0TE-01 

1582CE-03 

100006*01 

lOAC    CCKDITICN 

3TSPLAC6M6NTS 
JNT       >-CISP1. 

1  c.c 

2  -C.1C5E-C2 
2      -0. 1055-02 

*  O.C 
i         C.C 

icac  cooiTicik 

CTSPLACFHENTS 
J*T       X-CISPL 

1  C.C 

*  -0.226E-02 

2  -C. 2266-02 
+         C.C 

5         0.0 


1 

CEGPEE    OP   FREEDOM 

Y-CISPL 

O.C 

■0.275E-02 
■0.  1C5E-01 

O.C 

C.C 


2-01 SPl 

0.0 
-0.13  76-03 
-0.3*36-03 

0.0 

0.0 


DEGREE    CF    FREEDOM 


Y-CISFL 

O.C 
■0.271E-02 
■0. 161E-01 

C.C 

0.0 


2-DISPL 

0.0 
-0.  187E-02 
-0.7B*E-02 

0.0 

0.0 


*CT    AST    X 

0.0 

CO 

CO 

0.0 

0.0 


ROT   AST   X 
0.0 


0.0 
0.0 
0.0 

CO 


RCT    ABT   Y 

0.0 

0.206E-05 

0.2066-05 

0.0 

0.0 


ROT    ABT    Y 
0.0 

0.*16E-0* 

0.6886-0* 

0.0 

0.0 


ROT    A8T    Z 

0.0 
-0.5676-0* 
-0.8756-0* 

0.0 

0.0 


ROT   ABT    Z 

0.0 
-0.8*76-0* 
-0.1*3E-03 

0.0 

0.0 


T    ABT    I 

.0 

.  566*6E-0* 

.87*81E-0* 

.0 

.0 


T    ABT    Z 

.0 

.3*6896-0* 

.  1*26  <J  6-03 

.0 

.0 
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TABLE    IV 
EXAMPLE    1:       SAMPLE    CONSTRAINT    OUTPUT   INFORMATION 

DISPLACEMENT    CCNSTRAINTS    START    AT    NUMBER       2   ANO    STOP    AT    NUMBER       5 

TPUSS    ELEMENT    STRESS    CONSTRAINTS 

START    AT    NUMBER      6    ANO    STOP    AT    NUMBER    17 

FFAME    ELEMENT    STRESS    CONSTRAINTS 

START    AT    NUMBER    18    ANO    STOP    AT    NUMBER    41 

THE    NLJ^eER    CF    CCNSTFAINTS    NCTOT    =  41 

1)  -.4CESe2E-K2  -.161055E+O1  -.389445E+00  -.258895S+0C 

5)  -.174111E+C1  -.SS1588E+00  -.1018415+C1  -.943857E+00 

9)  -.1C1236E+C1  -.9?7t4CE+00  -.103769E+01  -. 102428E+C1 

12)  -.«75722E  +  CC  -.  102139E+Q1  -.388444E+00  -.111156E+01 

17  »  -,9C1716E-KC  -.  126755E+C1  -.732450E+CC  -.152875E+C1 

21)  -.471251E+C0  -.118559E+01  -.972695E+00  -.597988E+00 

25)  -.KC2C1E  +  CI  -.225C14E+CO  -.177499E+C1  -.127432E+01 

29)  -.95S641E+C0  -.  171444E  +  01  -.285557E+C0  -.lOOOOOE+Ol 

23)  -.9CS999E+C0  -.123411E+01  -.965558E+00  -.581741E-04 

37)  -.19«S994E  +  C1  -. S9SC9BE+G0  -. 10COOGE+ CI  -.132766E+01 

41)  -.951795E+CC 
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TABLE  V 
EXAMPLE  1:   FINAL  OPTIMIZATION  INFORMATION 


FINAL    OPTIMZATICN    INFCRMATION 

Cej   =         C.5E7  5C7E+C3 

OECISION    VARIABLES    (X-VECTORJ 

1)     C.348CJE+01       0.84992E+C1       0.48389E+01       0.21898E+00      0.15576E+00 

7)    C.45432E+02 

CCNSTPAINT    VALUES    (G-VECTORI 

l)-C.^03S£5+02  -0.161C6E+01    -0.  28945E  +  0C   -0. 25889E+00    -0. 17411E+01 

7)-C.1018«E+01  -0.94336E+00    -0.10124E+01    -0.98764E+00   -0.10377E+C1 

13)-C.975  72E+00  -0.10214E+01    -0.S8844E+CC    -0.111165+01    -0.901725  +  00 

19>-C.73245E+00  -0.  15 297E+C1    -0.47125E+CC    -0.11856E+C1    -0. 97270E+CC 

25 J-C.1402CE+01  -0.22501E+00    -0.17750E+01    -0  .12743 E+0 1    -0.95964E+QG 

21 >-G. 2855<E+O0  -C.  100C0E+C1    -C.  10000E+01    -0. 12341E+01    -O.C6556E+00 

27>-C.1999«E+01  -0.10000E+01    -0.10000E+01   -0.13277E+C1    -0.95180E+C0 

THERE   ARE  1    ACTIVE   CONSTRAINTS 

CONSTRAINT    NUMBERS    ARE 
it 

THEFE    ARE  0    VICLATED    CONSTRAINTS 

THERE   ARE  0    ACTIVE    SIDE    CONSTRAINTS 

TEFMNATICN  CFITEPION 

AeS(CBJ(  I  »-OBJ(  I-1JI     LESS    THAN    OABFUN    FOR    3    ITERATIONS 

NL*eEP    OF     ITERATICNS    =         12 

OBJECTIVE    FUNCTION    WAS    EVALUATED  114      TIMES 

CONSTRAINT    FUNCTIONS    *EPE    EVALUATED  114      TIMES 

THIS    RUN    PECUIREC       116      STRUCTURAL    ANALYSES 

NUKeER    OF    SECGNCS    REQUIRED    FOR    EXECUTIGN    =    2.79 

WEIGHT    OF    STRUCTURE    GIVEN    AREAS    S    LENGTHS 
WEIGHT=  0.5E751E+C3 

TCTAL    WEIGHT    INCLUDING    FIXED    MASSES 
TCTAL    WEIGHT*  C. 63751E+03 


JCINT 

CQCFOINATES 

JCINT 

X 

Y 

z 

1 

0.0 

0.0 

0.0 

2 

G. 1C0C0E+C3 

0.0 

0  .0 

2 

C.2CCC0E+C3 

0.0 

0.0 

4 

0.0 

0.84259E+02 

0  .0 

5 

O.C 

0.0 

0.45432E+02 

ELEMENT     INFORMATION    FOR    3AR    ELEMENTS 
ELEMENT-JCINT    RELATIONSHIPS 


ELEMENT 

NODE    1 

NOCE    2 

AREA 

LENGTH 

1 

2 

4 

C.3481E+C1 

0.  13C8E+03 

2 

2 

5 

0.8499E+01 

0.  1098E+03 

ELEMENT     INFORMATION    FOR    FRAME    ELEMENTS 
ELEfENT-JCINT    RELATIONSHIPS 

LNC    NCCE1    NCCE2  AREA  LENGTH  CHAR. DIM. 1       CHAR. DIM. 2 

LNC    NCCtl    ^"q^^E+OI       Q.1000E+03       0.*839E+01       0.2190E+00 

4  2  2      0l236EE+01       0.1000E+C3      0.4839E+01       0.1=>58E+00 
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V.   USER  GUIDE 

A.   INTRODUCTION 

In  developing  any  computer  code  for  engineering  analysis, 
it  is  necessary  to  additionally  develop  concise,  easily 
understood  documentation.   This  SADX  USER'S  GUIDE  is  written 
to  be  easily  understood  by  the  user  having  only  minimal  know- 
ledge of  the  FORTRAN  language .   The  format  follows  that  of 
the  optimization  code,  COPES/CONMIN  [Ref.  2]. 

This  chapter  is  devoted  to  acquainting  the  user  with  the 
code  and  necessary  input  data. 

3.   DESIGN  EXAMPLE 

The  simple  example  of  a  four-element  combined  truss /frame 
structure  is  used  to  demonstrate  some  of  the  features  of  the 
SADX  program. 

The  structure  is  shown  in  Figure  5.1,  and  consists  of 
two  tubular  frame  elements  along  the  x-axis  with  two  truss 
braces  to  The  y  and  z  axes  from  The  joint  between  the  frame 
elements.   A  non-structural  fixed  mass  is  attached  at  the 
outboard  end  of  the  second  frame  element  where  two  loads  are 
applied . 

C.   SADX  DATA 

The  SADX  program  reads  data  from  unit  5  and  writes  out- 
put on  unit  5.   Units  30  and  40  are  used  as  scratch  files. 
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ORIGINAL 


Figure  5.1   TRUSS-BRACED  CANTILEVER  BEAM 
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The  scratch  file  numbers  may  be  changed  by  changing  two  cards 
at  the  beginning  of  subroutine  SADX01.   The  SADX  program  has 
the  capability  to  read  unformatted  data.   The  following  sec- 
tions define  the  data  which  is  required  by  SADX.   The  data 
is  segmented  into  "BLOCKS"  for  convenience. 

SADX  data  begins  with  a  TITLE  card  and  ends  with  a  END 
card.   Comment  cards  may  be  inserted  anywhere  in  the  SADX 
data  stack  prior  to  the  END  card,  and  are  identified  by  a 
dollar  sign  ($)  in  column  1.   Data  may  be  of  either  the 
"110"  or  "F10.0"  type  or  may  be  free-format  separated  by 
commas  with  no  imbedded  blanks.   Lines  of  formatted  and 
unformatted  data  may  be  intermixed. 

1 .  Formatted  Data 

Formats  are  of  "110"  and  "F10.0"  type.   "I"  formats 
must  be  right  justified,  and  "F"  formats  must  have  the  deci- 
mal point.   The  number  of  cards  read  per  data  block  is  listed 
at  the  bottom  of  each  block. 

2 .  Unformatted  Data 

While  the  USER'S  MANUAL  data  sheet  defines  SADX  data 
in  formatted  fields  of  ten,  the  data  may  actually  be  read  in 
a  simplified  fashion  by  separating  data  by  commas  or  one  or 
more  blanks.   If  more  than  one  number  is  contained  on  an  un- 
formatted data  card,  a  comma  must  appear  somewhere  on  the  card. 
If  exponential  numbers  such  as  2.5+10  are  read  on  an  unfor- 
matted card,  there  must  be  no  imbedded  blanks.   Unformatted 
cards  may  be  intermingled  with  formatted  cards .   Real  numbers 
on  an  unformatted  card  must  have  a  decimal  point. 
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EXAMPLES 

Unformatted  data 

5,7,3.2 ,1.3  +  6,-5  ,2 

Equivalent  formatted  data 

col     10 20 30 40 5_0 50 70 80_ 

5       7     3.2   1.3+6      -5       0       2 

Unformatted  data 

2 

2,3 

2  3 

Equivalent  formatted  data 

col     10 20 30 40 50 60 7_0 30_ 

2 

2       3 

2  3 

NOTE:   The  third  line  of  data  contains  no  commas  and  is 
therefore  assumed  to  be  already  formatted. 

Placement  of  more  than  eight  unformatted  data  on  a 
card  will  create  two  (or  more)  formatted  cards  as  required 
Fields  of  zeros  will  be  created  if  more  data  are  required 
than  are  filled  on  an  unformatted  card. 

D.   CONSTRAINTS 

Constraints  are  calculated  and  stored  in  the  2    vector 

as  listed  in  the  following  chart.   The  total  number  of 

constraints 

NCON=NFREQ+2*MDSPLC+NLC*(2*NEB  +  4*NEF  +  NEB  +  2-NEF) 
freq     displ       stress    stress  buckl   buckl 
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Where  NCON  is  the  total  number  of  constraints,  NFREQ  is  the 
number  of  frequency  constraints,  NDSPLC  is  the  number  of 
displacement  constraints,  NLC  is  the  number  of  load  cases 
imposed,  NEB  is  the  number  of  bar  or  truss  elements,  and  MEF 
is  the  number  of  frame  elements .   When  any  of  the  constraints 
are  missing  from  the  G  vector,  all  constraints  are  moved  up 
in  the  vector.   For  example,  if  there  is  no  frequency  con- 
straint, then  a  displacement  constraint  would  fill  the  first 
location  of  the  G  vector. 

CONSTRAINTS  ARE  STORED  IN  THE  G  VECTOR  IN 
THE  FOLLOWING  ORDER: 

NFREQ  FREQUENCY  CONSTRAINTS 

2 -NDSPLC  JOINT  DISPLACEMENT  CONSTRAINTS 

STRESS  CONSTRAINTS  ARE  STORED  ELEMENT  BY  ELEMENT 

FOR  A  GIVEN  ELEMENT  CONSTRAINTS  ARE  STORED  BY  LOAD  CASE 

FOR  A  GIVEN  TRUSS  ELEMENT  AND  LOAD  CASE, 

CONSTRAINTS  ARE  STORED: 

TENSILE  STRESS  LOWER  LIMIT 

TENSILE  STRESS  UPPER  LIMIT 

EULER  BUCKLING  STRESS  LIMIT  (IF  APPLICABLE) 
FOR  A  GIVEN  FRAME  ELEMENT  AND  LOAD  CASE, 
CONSTRAINTS  ARE  STORED: 

NORMAL  STRESS  AT  LOW  NODE  LOWER  LIMIT 

NORMAL  STRESS  AT  LOW  NODE  UPPER  LIMIT 

NORMAL  STRESS  AT  HIGH  NODE  LOWER  LIMIT 
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NORMAL  STRESS  AT  HIGH  NODE  UPPER  LIMIT 
EULER  BUCKLING  STRESS  LIMIT  (IF  APPLICABLE) 
LOCAL  BUCKLING  STRESS  LIMIT  (IF  APPLICABLE) 

E .   EXAMPLE 

The  initial  configuration  of  the  braced  cantilever  beam 
is  shown  in  Figure  5.1  Stress  constraints  are  imposed  as 
well  as  constraints  on  Euler  and  local  buckling,  displace- 
ment, and  first  fundamental  frequency.   A  non-structural 
fixed  mass  is  applied  at  the  tip  of  the  beam,  and  two  load 
conditions  (PI,  and  P2 )  are  imposed.   The  structure's  own 
weight  will  be  considered  as  an  imposed  load  as  will  be  the 
fixed  mass . 

The  linking  of  design  variables  is  demonstrated  by  link- 
ing the  mean  diameters  of  the  two  frame  members  (Dl,  and 
D2  with  D1=D2).   Member  size  variables  are:   truss  member 
areas  (Al,  and  A2),  frame  member  mean  diameters  (Dl,  and 
D2 ) ,  and  frame  member  thicknesses  (Tl,  and  T2).   Member 
sizing  DESIGN  VARIABLES  are:   XA(i)=Al,  XA(2)=A2,  XA(3)=D1=D2 
XAU)=T1,  XA(5)=T2. 

Geometry  variables  are  the  attachment  points  (joints  4 
and  5)  of  the  two  truss  members  on  the  y  and  z  axes  (YU,Z5)  . 
Coordinate  DESIGN  VARIABLES  are:   XC(1)=Y4  and  XC(2)=Z5.   The 
structure's  weight  is  then  optimized  with  respect  to  member 
sizes  and  structure  geometry  for  a  total  of  seven  design 
variables . 
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1 •   Properties /Conditions 

Two  material  types  are  used:   type  1,  aluminum,  is 
used  for  the  truss  members;  and  type  2,  steel,  is  used  for 
the  frame  members.   The  weight  densities  of  the  materials 
(p)  are: 

type  1=0.1  lb/in3 

type  2=0.3  lb/in3 
The  Young's  moduli  of  the  materials  (E)  are: 

type  1   E  =  10.0E+5  psi 

type  2   E  =  29.CE+6  psi 
The  non-structural  fixed  mass  attached  at  the  tip  of  the 
structure  weights  250  lb.   The  applied  loads  are: 

PI      1000.0  lb  in  the  +y  direction 

P2      1000.0  lb  in  the  -z  direction 
The  acceptable  maximum  normal  stresses  are: 

ty-ce  1   -25000  psi  <  c     <  25000  osi 

J -  —  max  — 

tyDe  2   -36000  osi  <  a     <  35000  osi 

J  -  —   max  — 

Displacement  limits,  imposed  upon  joint  number  3  (the  tip) 
for  each  load  case  in  the  direction  of  loading,  are: 

load  case  1      y-direction      +/-  3.0  in. 

load  case  2      z-direction      +/-  3.5  in. 
3ounds ,  placed  on  the  positions  of  joints  4  and  5  along  the 
y  and  z  axes ,  are : 

joint  M-      y-coordinate      0.0  inches  to  200.0  inches 

joint  5      z-coordinate      0.0  inches  to  100.0  inches 
The  minimum  natural  frequency  of  the  structure  is  constrained 
to  be  greater  than  1Hz . 
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2 .   Input  Control  Parameters 

The  following  input  control  parameters  are  given 
for  ease  of  following  the  example: 


NEB  =  2 

NMT=2 
LBUCK=1 
NEIG1=2 
NSTRES=2 


NEF  =  2 
IDVCLC=2 

NFREQ=1 

NPRI=0 

NSTW=1 


NJ  =  5 
NDJ=2 
NFMASS=1 
NLC=2 

NFMW=1 


NCJ=3 

NEUBC=1 

NEIG=1 

NDSPLC=2 


Table  VI  is  a  listing  of  commonly  used  nomenclature . 

The  following  USER'S  MANUAL  is  divided  into  blocks 
A  through  P.   Appearing  directly  below  each  data  field  line 
are  the  parameters  for  the  TRUSS-BRACED  CANTILEVER  BEAM 
example .   It  is  important  to  note  that  the  user  may  choose 
any  consistent  system  of  units . 
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TABLE  VI 
COMMON  VARIABLE  NOMENCLATURE 


A  -    HEHBEB»S    CROSS -SECTIONAL    AREA 

BL  -    LO»ER    BOUND    ON    DISPLACEMENTS 

BO  -    UPPER    BOUND    ON    DISPLACEMENTS 

CHARDIM1, 2-    CHARACTERISTIC    DIMENSIONS    EOR    ERAMS    MEMBERS 

FOR    L5ECT.EQ.1:     MEAN    DIAMETER    AND    THICKNESS 
DIR  -    DIRECTION    1=X,    2=1,    3  =  Z,    4=rot    about    x, 

5=rot    about    y,    6=rot    about    z. 
S  -    YOUNGS    MODULUS 

EPSEIG  -    CONVERGENCE    TOLERANCE   OF    EIGENVALUE 

SOLUTIONS    f DEFAULT=.330  1) 
FC1,...FCN-    LOWER    BOUND    ON     FIRST, SECOND ,     ETC. 

SYSTEM    NATURAL    FREQUENCY 
FX  -    LOAD    FORCES    APPLIED    IN    THE    X    DIRECTION 

FY  -    LOAD    FORCES    APPLIED    IN    THE    Y    DIRECTION 

FZ  -    LOAD    FORCES    APPLIED    IN    THE    Z    DIRECTION 

GRAV  -    ACCELERATION    DUE    TO    GRAVITY     (DEFAULT=    386. tt) 

IDVCLC  -    THE    OPTIMIZATION    OPERATION    IDENTIFIER 

1  FOR    OPTIMUM    MEMBER    SIZE    ONLY 

2  FOR    BOTH    OPTIMUM    MEMBER    SIZE    AND    GEOMETRY 

3  FOR    OPTIMUM    GEOMETRY    ONLY 

IX  -    CONSTRAINT    IDENTIFIER.     IF    NON-ZERO    THE    X-DO? 

IS    CONSTRAINED 
IY  -    CONSTRAINT    IDENTIFIER.     IF    NON-ZERO    THE    Y-DCF 

IS    CONSTRAINED 
IZ  -    CONSTRAINT    IDENTIFIER.     IF    NON-ZERO    THE    Z-DCF 

IS    CONSTRAINED 
IXX  -    CONSTRAINT    IDENTIFIER.     IF    NON-ZERO    THE 

ROTATION    ABOUT    THE    X-AXIS    IS    CONSTRAINED 
IYY  -    CONSTRAINT    IDENTIFIER.     IF    NON-ZERO    THE 

ROTATION    ABOUT    THE    Y-AXIS    IS    CONSTRAINED 
IZZ  -    CONSTRAINT    IDENTIFIER.     IF    NON-ZERO    THE 

ROTATION    ABOUT    THE    Z-AXIS    IS    CONSTRAINED 
JN  -    JOINT    NUMBER    (GLOBAL) 

KEULER-    SULER    BUCKLING    COEFFICIENT    FOR    BAR 

ELEMENTS 
L3UCK  -    LOCAL    BUCKLING    CONSTRAINT    IDENTIFIER 

IF    LEOCK.NE.O,    LOCAL    BUCKLING    CONSTRAINTS 

SILL    BE    APPLIED    TO    TU3ULAR    FRAME     MEMBERS 
LC  -    LOAD    CONDITION 

LMASS  -    LUMPED    MASS    OPTIONS 

IF    LMASS.EQ.O    CONSISTENT    MASS    MATRIX    USED 

IF    LMASS. NE.O    LUMPED    MASS    MATRIX    USED 
LNO  -    ELEMENT    NUM3ER 

L5ECT  -    CROSS-SECTION    TYPE    IDENTIFIER 

L5ECT.EQ.1     INDICATES    A    TU3ULAR    MEMBER 
NCJ  -    THE    NUMBER    OF     CONSTRAINED    JOINTS 

NDJ  -    THE    NUMBER    OF    JOINTS    KITH     (OR    LINKED    TO) 

DESIGN    7ARIA3LES 
NDSG1  -    TRUSS    MEMBER    AREA    DESIGN    VARIABLE    NUMBER 

NDSG2  -    COORDINATE    DESIGN    VARIA3LE    NUMBER 

NDSG3  -    FRAME    MEMBER    CHARACTERISTIC    DIMENSION    1 

DESIGN    VARIABLE    NUMBER 
NDSGU  -    FRAME    MEMBER    CHARACTERISTIC    DIMENSION    2 

DESIGN    VARIABLE    NUM3ER 
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NDSPLC  -    THE    NUMBER   OF    DISPLACEMENT    CONSTRAINT    SETS 

NEB  -    THE    NUMBER    OF    TRUSS    ELEMENTS 

NSF  -    THE    NUMBER    OF     FRAME    ELEMENTS 

NEIG  -    NUMBER    OF    PRECISE    EIGENVALUES    TO    BE    FOUND 

NEIG1  -    NUMEER    OF    SI3E NVALUE/EIGEN VECTOR    TO    BE 

EVALUATED    DEFAULT    =MIN(2*NEIG,    NSIG+8) 
NEUBC  -    EULER    BUCKLING    CONSTRAINT    IDENTIFIERS. 

IF    NEUBC. MS. 0    EULER    BUCKLING    CONSTRAINTS 

WILL    EE    IMPOSED    ON    BAR    AND    FRAME    ELEMENTS 
NFMASS  -    NUMBER   OF    NON- STRUCTURAL    FIXED    MASSES 

ATTACHED    TO    THE    STRUCTURE 
NFMW  -    FIXED    MASS    WEIGHT    IDENTIFIER 

IFfNFMW.NE. 0)     FIXED    MASSES    WILL    BE 

CONSTDERSD    AS    LOADS 
NFREQ  -    NUMEER    OF    FREQUENCY    CONSTRAINTS 

NID  -    NUMBER    OF    INDEPENDENT    DEGRESS    OF    FREEDOM 

NJ  -    THE    NUMBER    OF    JOINTS 

NLC  -    NUMBER    OF    LOADING   CONDTIONS 

NLJ  -    NUMBER    OF    LOADED    JOINTS    FOR    THIS    LOAD    CONDITION 

NMT  -    NUMBER    OF   SEPARATE    MATERIAL    TYPES 

NPRI  -    INPUT    PRINT    CONTROL 

IF    NPRI.N3.0    INPUT    VALUES    are    NOT    PRINTED 

IF    NPRI.SQ.5    LOCI/LOCR/IA/RA    ARRAYS    PRINTED 
NSTRES  -    THE    FORCE/STRESS    PRINT    IDENTIFIER 

I?    NSTRES. SQ.O    MEMBER    FORCES    PRINTED 

IF    NSTRES. SQ.1    MEMBER    STRESSES    PRINTED 

IF    NSTRES. EQ. 2    BOTH    are    PRINTED 
NSTW  -    STRUCTURE   WSI3HT    IDENTIFIER 

iF(NSTW. NE. 0)     THE    STRUCTURE'S    OWN    WEIGHT 

WILL    BE    CONSIDERED    AS    LOADS 
POISSN  -    POISSON'S    RATIO 

PHO  -    MATERIAL    WEIGHT    DENSITY 

SIGMIN  -    MINIMUM    ALLOWABLE    NORMAL    STRESS 

SIGMAX  -    MAXIMUM    ALLOWABLE    NORMAL     STRESS 

TX  -    TORSIONAL    MOMENT    APPLIED    ABOUT    THE    X-AXIS 

TY  -    TORSIONAL    MOMENT    APPLIED    ABOUT    THE    Y-AXIS 

TZ  -    TORSIONAL    MOMENT    APPLIED    ABOUT    THE    Z-AXIS 

XA  -    INITIAL    VALUE     OF    AREA    DESIGN    VARIABLE 

XAL  -    LOWER    BOUNDS    0 N    XA 

XAU  -    UPPER    BOUNDS    ON    XA 

XC  -    INITIAL    VALUE    C?    COORDINATE    DESIGN    VARIA3LE 

XCL  -    LOWER    BOUNDS    0  N    XC 

XCU  -    UPPER    BOUNDS    0  N    XC 


i+O 


DATA    BLOCK       A 


DESCRIPTION:    Title    Card 


Format   and  Example 


:itlz 


FORMAT 
20A4 


TRUSS-BRACED    CANTILEVER    BEAM 


FIELD 


CONTENTS 


1  ANY    80    CHARACTER    TITLE    MAY    BE    GIVEN    ON    THIS 

LINE 
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DATA    BLOCK      B 


DESCRIPTION:    Control   Parameters 


Format   and   Example 


NEB 

NEF 

NJ 

NCJ 

HHI 

IDVCLC 

NDJ 

FORMAT 
7110 

NEUBC    LBUCK    NFREC    NFMASS     MEIG       NSTG1     NPRI 


orma- 
7110 


1 

1 

1 

1 

1 

0 

NLC 

NDSPLC 

NSTHES 

Hsrw 

NF  M  H 

FORMAT 
5110 

NOTE:  DEFINITIONS  OF  PROGRAM  INPUT  CONTROL 
PARAMETERS  ARE  LISTED  ON  NEXT  PAGE 
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FIELD  CONTENT 

FIRST  CARD 

1  NEB-number  of  bar  elements 

2  NEF-number  of  frame    elements 

3  NJ-number   of   joints 

4  NCJ-number   of  constrained    joints 

5  NHT-number  of  seperate  material   types 

6  IDVCLC-design  variable   control   parameter 

If     (IDVCLC.EQ.1)     NDV=NDVAR1 
optimizes    wrt    member   size    only 
If     (IDVCLC.EQ.2)     NDV=NDVAR1+    NDVAR2 
optimizes   wrt    member   size   and   geometry 
If     (IDVCLC.EQ.3)     NDV=NDVAR2 
optimizes    wrt    geometry   only 
NDJ-total    linked    and    design    variable    joints 
(i.e.    number    cf    'movable1    joints) 

SECOND    CARD 

1  NEUBC-Euler    buckling   constrain-    ideatifier 

If    (N20EC.NS.  0)    EULER    buckling    constraints 
will    be    imposed  on    bar    elements 

2  LBUCK-local   buckling   constraint   identifier 

If    (LBUCK.NE.O)    local    buckling    constraints 
will    be    imposed   on    tubular    frame    elements 

3  NFREQ-number    of    freguency    constraints 

4  NFMASS-number   of   fixed    masses   attached   to 

structure 

5  NEIG-number   of    precise   eigenvalues   to 

be   evaluated 
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6  NEIG1-number   of    eigenvalues   to    be    evaluated 

DS?ADLT=min.     of    (2*NEIG    ,       NEIG  +  3) 

7  NPR1-prirt  control  identifier 

NPRI.ne.O  input  info  not  printed 
NPR1.eq.5  RA/IA/LOCR/LDCI  arrays 
will    be    printed  for   debugging 

THIRD    CARD 

1  NLC-number   of  load  conditions 

2  NDSPLC- number   of   displacement   constraints 

3  NSTRES-fcrce/st ress    print    identifier 

If    (HSTRES.EQ .0)     stresses   are   printed 
If    (NSTRES.SQ . 1)    forces    are    printed 
If    (NSTRES.SQ .2)    both    ire    printed 

4  NSTW-structure    weight   identifier 

If    (NSTW.NS.O)    rhe    structure's    weight    is 
considered   is    loads 

5  NFMW-fixed  mass    weight    identifier 

If    (HFMW.NE.O)    the    fixsd    masses   are 
considered   as    loads 
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DATA    BLOCK      C 


DESCRIPTION:    Dynamic    Analysis    Information 


Format   and  Example 


LMASS 

GRAV 

EPSEIG 

FORHAI 

I10,2?10.0 

"1 — 

386.4 

0.0 

FIELD  CONTENTS 

1  LMASS- lumped    mass   option     (if    LMASS. 52.0) 

the   lumped    mass    matrix    is    used. 

2  GRAV-accleration   due    to    gravity 

(default=386.  4    inches/sec2) 

3  EPSEIG-convargence   tolerance    on   eigenvalue 

solution.     (dafault=.0031) 
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DATA    BLOCK       D 


DESCRIPTION:      Joint   Coordinates 


Format   and   Example 


JN 

X 

Y 

Z 

FORMAT 
110, 3F10.0 

1 

0.0 

0.  0 

0.0 

2 

100.0 

0.0 

0.0 

3 

2  00.0 

0.0 

0.0 

t) 

0.  0 

150.0 

0.0 

5 

0.0 

0.  0 

50.0 

FIELD 


CONTENTS 


1  JN- joint    cooriina~e    number 

2  X-x    coordinate 

3  Y-y    coordinate 
ix  Z-z    coord  inate 


NOTE:    Number   of    cards   reai=NJ 
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DAT!    BLOCK       5 

Omit    this    block    if   NDJ=0    in    BLOCK    B 
DESCRIPTION:      Coordinate    Design   Variable   Linking    Data 
Format   and   Example 


JN 

IX 

II 

IZ 

PCX 

PCI 

PCZ 

FORMAT 
4i10,3f 10.0 

4 

0 

1 

0 

1.0 

1.0 

1  .0 

5 

0 

0 

2 

1.0 

1.0 

1.0 

FIELD  CONTENTS 

1  IX-desigr    variable   associated  with   x    coori 

2  lY-desigr.   variable   associated  with  y    coord, 

3  IZ-desigr.    variable   associated  with   z    cccr3 

4  PCX-participation   coefficient  of   x-ccord. 

5  PCY-participati en   coefficient  of   y-coord. 

6  PCZ-Darticioati on   coefficient  cf   z-cocrd. 


NOTE:    Number   of    cards   read=>TDJ 
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DATA    BLOCK      F 

21 SCRIPT ION:    Material   Properties 

Format   and   Example 


E 

RHO 

SIGMIN 

SIGMAX 

KEULER 

POISSN 

FORMAT 
6F10.0 

1. OE+7 

0.1 

-25000. 

25000. 

2.3 

0.27 

2.9E+7 

0.3 

-360  00. 

3  60  00. 

2.0 

0.27 

FIELD  CONTENTS 

1  E-Young's    Modulus 

2  RHO-marerial    density 

3  SIGMIN- minim  urn    allowable   stress 
a  SIGMAX-maximum    allowable    stress 

5  KEULSR-Euler   buckling    coefficien- 

6  POISSN-Pcisson1 s   ratio 


NOTE:    Number  of   cards    rsad=NMT 


(+8 


DATA    BLOCK      G 


Omit    this    block    if    NE3=0    in    BLOCK    3 


DESCRIPTION :    Bar   Element   Information 


Format   and   Example 


LNC 


NODE1 


NODE2 


MAI COD 


NDSG1 


LSECT 


FORMAT 

5110 
f 10,110 


3 

2 

4 

-1 
1 

1 

3.0 

1 

2 

2 

5 

1 

2 

3.0 

1 

!LD 


CONTENTS 


1  LNO-element   number 

2  NODEl-glcbai   number   associated    with   the 

element's    1st    node 

3  NODE2-giobal   number   associated    with   the 

element's   2nl    node 

4  :*  AT  COD- mater  ial    type    of   this    element 

5  HDSG1 -design    variable    number    associated    with 
this    element's   area 

6  A-m ember   cross- sect ional    area 

7  LSECT-cross-sscticn   type    identifier 

LSECT. EQ.1    indicates    a    tubular   member 

NOTE:    Number   of  cards   rsad=SEB 
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DATA    3L0CK       H 

Omit    this    block    if    N3F  =  0    in    BLOCK    B 
DESCRIPTION:         Frame   Element      Information   Format 
Example 


and 


INC 


NODE1 


NODE2 


MAT COD 


NDSG3 


NDSG4 


LSECT 


FORMAT 
7110 


3 

1 

2 

2 

2 

4 

1 

CHARDIM1 

CHARDIM2 

FORMAT 
2F10 

5.0 


1.0 


u 

2 

3 

2 

3 

5 

1 

5.0 

1.0 

;IELD  CONTENTS 

1  LNO-element   number 

2  NODS1-global   number   associated    with    the 

element's    1st    node 

3  NODE2-global  number   associated    with   the 

element's    2nd    node 

4  MATCOD-material    type    of    this    element 

5  NDSG3-design    variable   number   associated    with 

the   element's    1st    characteristic   dimension 

6  NDSGU-design    variable   number   associated    with 

the    element's    2nd   characteristic   dimension 

7  LSECT-crcss-ssction   type    identifier 

LSECT. SQ.1    indicates    a    tubular    member 
NOTE:    Number    of   cards   read=NSF 
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DATA    BLOCK      I 


DSSC3IPTICN :    Joint   Constraint   Information 


Format   and   Example 


JN 

IX 

II 

IZ 

IXX 

IYY 

IZZ 

FORH&T 

7110 

1 

1 

1 

1 

1 

1 

1 

<L 

0 

0 

0 

0 

0 

0 

3 

0 

0 

3 

0 

0 

0 

n 

1 

1 

1 

1 

1 

1 

5 

1 

1 

1 

1 

1 

1 

FIELD 


CONTENTS 


1 
2 
3 

a 
5 

6 

7 


J H-  joint    number 

IX-    x-dis  placement   constraint    identifie: 

II-    v-dis placement    constraint    in- 

IZ-    z-displacan  ent   constrain4- 


II    -      _u< 


IXX-    x-axis   rotation   constraint   identirrer 
III-   y-axis   rotation  constraint    identifier 
IZZ-   z-axis   rotation  constraint   identifier 
if   . NE. 0   -   corresponding    DO?    constrained 


NOTE:    Number   of  car  is   rsad-NCJ 
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DATA    2L0CK       J 

Omit      this        block      if      HLC-0 
OISCPIPTION:    Joint   Loading   Information 

Format   and  Example 


in         3  LOCK      3, 


NLJ 

FOPilAT 

no 

JN 

FX 

FY 

F2 

TX 

IY 

TZ 

POHH&T 

I10,6F10 

2 

0.0 

1C00. 

0.3 

0.0 

0-0 

G.  0 

2 

0.0 

0.0 

-1000  . 

0.0 

3.  0 

0.0 

FIELD 
1 


CONTENT 


NLJ-number  of   leaded    joints    for   this    load 
condition 

1  JN- joint    number 

2  FX 

3  FY-   Forces  in  the  X,Y,Z    iirections 

4  FZ 

5  TX 

6  TY-    Moments  about    the   X,Y,Z    axes 

7  TZ 

NOT2:    Number   of  cards   read   per    sex=NLJ 
Number   of  sets    of  cards   read=NLC 
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DATA    3 LOCK       K 


Omit    this    block    if    NFMASS=0    in    BLOCK    B 


DESCRIPTION:    Fixed    Mass    Information 


Format   and   Example 


JN 

MASS 

FORMAT 
110, F10 

250.0 


FIELD 


CONTENTS 


1  JN- joint    number 

2  HASS-ooint   mass    at    joint    (JN)     in    force   units 


NOTE:    Number   of  cards   read=NFMASS 
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DATA    BLOCK      L 

Omit    this    block    if  IDVCLC=3 

DESCRIPTION:    Design  Variable    Information 

(MEMBER    SIZE    VARIABLES) 

Format    and  Example 

XA(1) 

XA(2) 

....     XA(NDVAR1) 

FORMAT 
8F10.0 

20.0 

20.0 

20.0 

2.0 

2.0 

XAL  (1) 

XAL  (2) 

.     .     .     .  XAL(NDVARI) 

FORMAT 
8F10.0 

0.50 

0.50 

4.0 

0.10 

0.  10 

XAU  (1  ) 

XAO  (2) 

.     .     .     .XAU(NDVARI) 

FORMAT 
8F10.0 

25.0              35.0 

25.0               2.5 

4.0 

FIELD                                                     CONTENTS 

1           XA- initial  value   of    area    design    variables 

2          XAL-lower   bounds   on    area    design   variables 

3          XAU-upper   bounds   on    area    design   variables 

NOTE:    read  one   value   of    XAr    XAL,    XAU   for    each 

independent    member    size    variable    defined    in 

BLOCKS    G    and    H 

Number    of   cards    read   =    as    required 
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DATA    BLOCK       H 

Omit    this    block    if   IDVCLC=1 

DESCRIPTION:    Design  Variable    Information 

(COORDINATE    VARIABLES) 

Format   and   Example 


XC  (1) 

XC(2) 

....     XC(NDVAR2) 

FORMAT 
8F10.3 

150.0 


50.0 


XCL  (1) 

XCL(2) 

.     .     .      .XCL(NDVA32) 

FORMAT 
3F10.0 

0.0 

0.0 

XCD  (1) 

XCU(2) 

FORMAT 
8F10.0 

2  00.0 


100.0 


FIFLD 

1 

2 
3 


CONTENTS 

XC-initial  value  cf  ccord.  design  variables 
XCL-lower  bounds  on  coord,  design  variables 
XCU-upper   bounds   on    coord,    design   variables 

NOTE:    read   one   value   of   XC,XCL,XCU    for    each 
independent    coord,    variable    defined   in    3L0CK    D 

Number    of   cards    read    =as    required. 
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DATA    BLOCK      N 


Omit    this    block    if    NDSPLC=0    in    BLOCK    B 
DESCRIPTION:    Joint    Displacement   Constraint   Information 


Format   and   Example 


JN 

DIB 

LC 

3L 

30 

FOFttAT 
3I10,2F10.0 

3 

2 

1 

-3.3 

3.0 

3 

3 

2 

-3.5 

3.5 

FIELD 


CONTENTS 


3 
4 

5 


JN- joint    numbsr 
DIH-direction 

1=x,2=y,3=z   displacement 

4=x,5=y,S=z   rotation 
LC-load   condition 
3L-Iower    bound    on  displac^m^n- 
BJ-upper  bound    on  displacement 


NOTE:    Number   of  cards   read=   NDSPLC 
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DATA    BLOCK      0 

Oait    this    block    if    HPRBQ=0    in    3LOCK    3 
DESCRIPTION:    Frequency   Constrain-    Information 
Format   and   Example 


FC1 

FC2 

?C3 

FOHHAT 
8F10.0 

1.0 


CONTENTS 

FC1 -   lower  bound  on   first   natural   frequency 

constraint   in    Hz.     {cycles    per    second) 
FCN-   lower  bound   en    NF3E2-th   natural 
frequency   constraint.    Id    Hz. 
(cycles    per    second) 


NOTE:    Number   of  cards   read   =    is    require! 
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DATA    BLOCK       P 


DESCRIPTION:    End   card 


Format   and   Example 


END 

FORMAT 
3A1 

END 


NOTE:  This  card  BUST  appear  at  the  end  of  the 
SADX  data. 
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VI.   NUMERICAL  EXAMPLES 

A.   INTRODUCTION 

Design  of  three-dimensional  truss  and  frame  structures 
are  presented  herein  and  the  corresponding  numerical  results 
are  summarized  to  demonstrate  the  use  of  the  SADX  code. 

The  examples  begin  with  the  truss-braced  cantiliver 
beam. 

3.   EXAMPLE  1:   TRUSS-3RACED  CANTILIVER  BEAM 

The  simple  truss-braced  cantiliver  beam,  as  shown  in 
Figure  6.1,  has  been  previously  used  for  the  SADX  USER'S 
MANUAL  example.   It  consists  of  two  steel  tubular  frame 
members  with  a  common  outer  diameter  and  different  wall 
thicknesses  arranged  as  a  cantilever  beam  along  the  X-axis. 
There  is  a  fixed  250  pound  mass  at  the  tip  of  the  beam. 
Two  aluminum  truss  members  are  attached  from  the  beam  mid- 
point to  points  en  the  Y  and  Z  axes .   This  structure  is 
designed  for  optimum  member  size  and  geometry  under  a  set 
of  two  load  conditions  and  subject  to  constraints  on  maximum 
stress,  maximum  joint  displacement,  Euler  and  local  buckling, 
and  minimum  structure  natural  frequencies .   The  weight  of 
the  non-structural  fixed  mass  and  the  structure's  own  weight 
are  imposed  as  loads.   There  are  five  member  size  and  two 
coordinate  design  variables,  and  a  total  of  41  constraints. 
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ORIGINAL 


Figure  5.1   TRUSS-BRACED  CANTILEVER  BEAM 


6  : 


The  number  of  analyses  required  for  this  design  is  116  using 
2.75  seconds  of  CPU  time  and  terminating  on  the  12th  itera- 
tion.  Of  the  analyses  conducted,  84  were  required  for  the 
calculation  of  gradients.   The  weight  of  the  structure, 
including  non-structural  fixed  masses,  is  reduced  from  9542 
pounds  to  333  pounds.   Results  are  given  in  Table  VII. 

C.   EXAMPLE  2:   TWO-TIER  3-D  PORTAL  FRAME  WITH  TRUSS  X-BRACES 

The  two-tier  three-dimensional  portal  frame  with  truss 
member  diagonal  braces,  as  shown  in  Figure  6.2,  is  a  symme- 
tric moderately  sized  case  that  can  be  analyzed  easily  by  a 
variety  of  other  codes.   There  are  20  truss  elements,  15 
frame  elements ,  and  four  non-structural  fixed  masses .   The 
material  used  is  steel.   This  structure  is  designed  for 
optimum  member  size  and  geometry  under  a  set  of  three  lead 
conditions  and  subject  to  constraints  on  maximum  stress, 
maximum  join-:  displacement,  Euler  and  local  buckling,  and 
minimum  structure  natural  frequencies .   The  weight  zz    the 
non-structural  fixed  masses  and  the  structure's  own  weight 
are  imoosed  as  loads .   There  are  10  member  size  and  five 
coordinate  design  variables  and  a  total  of  493  constraints. 
The  number  of  analyses  required  for  this  lesign  is  2  35  using 
120  seconds  of  CPU  time  and  terminating  on  the  21st  itera- 
tion.  Of  the  analyses  conducted,  30  5  were  required  for  the 
calculation  of  gradients.   The  weight  of  the  structure, 
including  non- structural  fixed  masses,  is  reduced  from  93CZ 
pounds  to  1462  pounds.   Results  are  given  in  Table  VIII. 
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EXAMPLE  1:   FINAL  3PTIMIZATI0N  INFORMATION 

FINU    OPTIMIZATION    INFORMATION 

G8w    =         C.5875CTE  +  C3 

THEFE    ARE  1    ACTIVE    CONSTRAINTS 

CONSTRAINT    NUMBEFS    APE 
26 

THEPE    ARE  0    VIOLATED    CONSTRAINTS 

Tl-EPE    ARE  C    ACTIVE    SIDE    CONSTRAINTS 

TSWINATICN    C^ITEF  ICN 

A£S(CBJ  (  I  l-CBJ(  1-1  ))     LESS    THAN    CABFUN    FOR    2    ITERATIONS 

NLMEEP    OF    ITERATIONS    =         12 

OBJECTIVE    FUNCTION    WAS    EVALUATED  114      TIMES 

CONSTRAINT    FUNCTIONS    fcERE    5VALUATEO  114       TIMES 

THIS    PUN    FECUIRED       116      STPUCT'JPAL    ANALYSES 

NUMBER    OF    SECONCS    REQUIRED    FOR     EXECUTION    =    2.79 

WEIGHT    OF    STRLCTUPE    GIVEN    AREAS    5    LENGTHS 
WEIGHT  =  G.5S751E+C3 

TOTAL    -EIGHT    IftCLLDING    FIXED    MASSES 
TOTAL    WEICHT=  C83751E+03 

JOINT    COCFCINATSS 
JOINT  X 

1  0    c 

2  C1CCCOS-K3 

3  G.2CCC3E+C3 

4  0.0 

5  0.0 

ELEMENT    INFORMATION    FOR    2AR    ELEMENTS 

ELECENT-JCINT    RELATIONSHIPS 

ELEMENT      NODE    1       NODE    2  AREA  LENGTH 

1  2  *  :.2^81E  +  C1  3.  12C3E+03 

2  2  5  D.8499E+01  3.1Q53E+03 

ELEMENT     INFORMATION    FOR    FRAME    ELEMENTS 
ELEVENT-JC!NT    RELATIONSHIPS 

LKC    KC0E1    NCCE2  *RE*  LENGTH  Ct-AR.DM.l      CHAR. DIM. 2 

-5  12       3.222<S+01       3.i3G0c+C2       3.*82Qr>01       Q.2190E+00 

4  2  2      0.226SE+C1      0.1GCGE+C3      0 . 4g25E+Gl       3.155BE  +  GO 

EIGENVALUES    AND    EIGENVECTORS 

NCDE    NUMBER  1 

FREQUEKCY      =     C.41S8EE+C2    CPS 

C.2621<E+C2    RPS 
EIGENVALUE    =    C.c92"'CE  +  C5 
EIGENVECTOR       DEGREE    CF    FREEDOM 
JOINT  1  2  3 

1  CO  CO  0.0 

2  CO  C.  C  CO 

~  -G.2342SE-C1      -C32985E-01       -C.7^5C7E-G1 

4  -0.50777E-02      -C.10745E-01         0.15320E-03 

5  -C22A25E-CI      -C.  17155E-01         C130CCE+01 


Y 

I 

0.0 

0  .0 

:.  c 

O.C 

0.0 

0.0 

C  2^259E+02 

3  .0 

CO 

3.-5432E+02 

:  : 


EXAMPLE  2 


Figure  5.2   TWO-TIER  3-D  PORTAL  FRAME  WITH  TRUSS  X-3RACF5 
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7.1.511  VIII 

example  2:  final  optimization  information 

FIML    OPTIMIZATION    INFORMATION 

Cej    =  C.144ie2E+C4 

TFEPE    AP  6  3    iCTIVE    CONSTRAINTS 

CONSTRAINT    NUMBERS    ARE 
46  91       196 

THERE    ARE  C    VIOLATED    CONSTRAINTS 

ThEFF    ARE  1    ACTIVE    SIDE    CONSTRAINTS 

DECISION    \APIABLE5    AT    LOWER    OS    LPPER     BCtMOS 
CINUS    INCICATES    LOWER    SOUND) 
-14 

TERMINATION    CRITERION 

AES(C5J(  I)-CBJ(I-U)    LESS    THAN    CAB  FUN    FOR    3    ITERATIONS 

NUMBER    OF    ITERATIONS    =         21 

OBJECTIVE    FLNCTICN    WAS    EVALUATEC  383      TIMES 

CONSTRAINT    FUNCTIONS    WERE    EVALUATED  333       TIMES 

TUS    PUN    PECUIREC      385       STRUCTURAL    ANALYSES 

NLPeER    OF    SECCNCS    FECLIREC    FOP     EXECUTION    =    120.7 

WEIGHT    CF    STRUCTURE    GIVEN    AREAS     5    LENGTHS 
WEIGHT=  0.1441cE+C4 

:CING    FIXED    MASSES 
C. 1641SE+04 

Y  Z 

3.0  0.9C0C2E+02 

CO  -3.9G3C2E*02 

CO  -0.9C3C2E+02 
CO  3.9C3C2E  +  02 

C20C30E+02  0.69074E+02 

CZ0GCGE  +  C2  -3.6937^5+02 

:.2C000E+O2  -0  .69G7-»E+02 
C.200Q0E+02  3.69074E+02 

C  100  00  E +03  3.120C0E  +  O3 

3.IO030S+O2  -0.1GOOOE+03 

C.  100C0E  +  03  -J.i:0C3E+03 
C10000E+C3  0.1G3CCS+03 

ile^ent   information  fcr   bar  elements 
:  ie  pe*t- joint  felaticnships 

:lement     node  i     ncce  2  area 

116  3.<+255E  +  C0 

2  2  3.4255E+C0 

3  2  7  G.4255E+CC 

4  3  6  0.4255E+CG 

5  2  8  0.4255E+CO 

6  4  7  0.4255E+C0 

7  4  5  0.H255E+C0  3.1^c5E+03 

8  18  C4Z55E+CC  3.  1465E+03 

9  5  7  0.4255E+C0  G.2U3E+03 

10  6  8  0.4255E+00  0.2113E+O3 

11  5  1C  3.  154^S*Ci  0.  1381E  +  03 

12  6  9  C154-E+01  0.1381E+03 
12  6  11  0.154hE+C1  3.1996S+03 


-:- 


TCT 

AL    WEIGHT    INCL 

TCTAL    *EIGHT= 

JCI 

hi   CGCFDINATES 

INT 

X 

1 

0.62385E+C2 

2 

0. 62285E+C2 

■a 

-0. 622S5E+C2 

A 

-0.62235  E  +  C2 

5 

0. 8C229E+C2 

3.8C239E+C2 

7 

-0.8C239  E+02 

8 

-0. SC229E+C2 

9 

0. 1C2C0E  +  C3 

1C 

0.1CCC0E+C2 

i: 

-0.  1C00OE+C3 

12 

-0.  1CCC0E+C3 

D.   EXAMPLE  3:   DD-9  6  3  FOREMAST 

Example  three  is  a  redesign  of  the  forward  mast  on  the 
DD-9  6  3  of  SPRUAMCE  class  destroyer.   The  DD-9  6  3  foremast  has 
been  chosen  as  the  third  test  case  for  the  following  reasons : 
1)  the  structure  is  typical  of  the  masts  found  on  many  com- 
batants in  the  United  States  and  other  navies,  2)  high  top- 
side weight  reduction  is  desirable  from  a  stability  viewpoint 
for  any  ship,  3)  the  structural  members  are  predominantly 
tubular,  4)  the  problem  is  sufficiently  complex  to  make 
conventional  design  methods  cumbersome,  5)  the  member  size 
and  loading  information  is  available. 

The  structure  as  shown  in  Figures  6.3  through  5.5  is 
constructed  of  172  aluminum  frame  members .   The  material 
used  in  the  structure  is  50  36  H32  aluminum,  an  alloy  with 
moderate  strength,  good  weldability,  and  good  corrosion 
resistance . 

This  structure  is  designed  for  optimum  member  size  under 
a  single  load  condition  and  subject  to  constraints  on  maxi- 
mum member  stress,  maximum  joint  displacement,  Zuler  and 
local  buckling. 

Some  structural  simplifications  are  made.   The  weights 
of  mast-mounted  radars,  antennas,  and  platforms  are  modeled 
by  a  series  of  fixed  masses  which  are  imposed  as  loads  along 
with  the  structure's  own  weight.   The  forces  due  to  halyards 
and  wire  antennas  are  applied  as  loads .   Inertial  forces  due 
to  ships  motion  are  calculated  for  the  initial  design  point 
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Figure  5.3   DD-363  FOREMAST  PORT  SIDE 
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Figure  5.4   DD-963  FOREMAST  STARBOARD  3i:i 
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Figure  5.5   DD-3S-3  FOREMAST  FORWARD  SIDE 
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Figure    5.6       DD-963    FOREMAST   AFT    £-11 
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con  sis":    of    The    product    of    a   member's 
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ccorcir.a  _e -dependent  load— f  a crer  eerliei  at  each  end  of  th< 
element  [Ref.  6].   The  load-factor  multiplier  is  dependent 
on  the  distances  from  a  ecir. t  e~  The  i — tersection  z~    the 
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TABLE    IX 

EXAMPLE    3:       FINAL   OPTIMIZATION    INFORMATION 

FIML    OPTIMIZATION    INFORMATION 

CEJ    =  C.3986C4E+05 

JHEPE    ARE  2    ACTIVE    CONSTRAINTS 

CCNSTRAINT   NUMBERS    ARE 
660       866 

TKEPE    ARE  0    VICLATEO    CONSTRAINTS 

THERE    ARE  C    ACTIVE    SIDE    CONSTRAINTS 

TEFMINATICN    CPITEPICN 

ABS(OBJ(  I)-08J(I-1)     LESS    THAN    CABFUN    FOR    3    ITERATIONS 

NUMBER    CF     ITERATIONS    =         17 

OBJECTIVE    FUNCTION    WAS    EVALUATED  553  TIMES 

CCNSTRAIM    FUNCTIONS    WERE    EVALUATSO  553  TIMES 

THIS    RUN    FEQUIREC       555       STRUCTURAL   ANALYSES 

WEIGHT    OF    STRUCTLRE    GIVEN    AREAS    £    LENGTHS 
WEIGHT=  0.29860E+05 

TOTAL    WEIGHT     INCLUDING    FIXED    MASSES 
TCTAL    WEIGHT=  0.36746E+05 

NUMBER    CF    SECONDS    REQUIRED    FOR    EXECUTION    =   4482.46 

hC.    CF    FFAME    ELEMENTS    =      172 
ELEMENT-JOINT    RELATIONSHIPS 

LNC    NCCE1    N0CE2            AREA                    LENGTH  CHAR. DIM. 1  CHAR.0IM.2 

1  4            12       C.5S7SE  +  02       0.1061E+03  C.19C05+02  0.1002E+01 

2  3             11      0.597SE+C2      0.1061E+03  C.19C0E+02  C.1CC2E+01 
2                            S       C.5979E+C2       0.H12E+03  C.1900E+02  0.1002E+01 

4  2            10       C.597SE+C2       Q.1112E+03  C. 1900E+02  0.10C2E+01 

5  4               8      0.6247E+01       0.1342E+03  C.6562E+01  0.302CE+00 

6  3               6       C.6247E  +  C1       0.1342E+03  C.6563E+C1  0.3O30E+OO 

7  18       0.6247E+01       Q.1369E+C3  C.6563E+01  C.2C3CE+00 

8  2               6       0.6247E+01       0.136<=>E-K>3  0.6563E+01  0.3030E+00 

9  8             12       C.62475+C1       0.1295E+03  C.6562E+01  0.3030E+00 

10  6             11       0.6247E+01       0.1295E+03  C.o563E+01  0.202CE+OO 

11  6               9       C.6247E+G1       0.11326+03  G.6562S+01  0.3030E+00 

12  6            10      0.6247E+C1       C.1132E+03  C.6563E+01  C.3C3CE+00 

13  3               7       0.6247E+01       3.6846E+02  0.6563E+01  0.303CE+00 
1«            4               7       C.6247E+C1       0.6846E+02  C.6563E+01  0.3030E+00 

15  7             11       0.6247E+C1       0.6346E+02  C.6563E  +  C1  C.3C3GE+00 

16  7             12       0.62<WE+01       0.6846E+02  G.6563E+01  0.3030E+00 

17  1                5       C.6247E+G1       0.I196E+03  C.6563E+01  0.3030E+00 

18  2               5       0.62^7E+C1       0.1196E+03  C.6563E+0I  0.3G2CE+00 
IS            5               9       C.6247E+C1       0.1Q97E+03  0.6563E+01  0.3020E+00 

20  5             1C       C.6247E+C1       Q.1097E+03  C.6563E+G1  0.3C3CE+00 

21  9            12      Q.6247E+Q1       0.2172E+03  Q.6563E+01  0.2C205+00 

22  1C             11       C.6247E  +  C1       0.2172E+03  C.6563E+01  0.3030E+00 

23  11             12       0.6247E+G1       G.84C0E+02  C. 6563E+Q1  G.3C3CE+00 

24  c            1C       0.6247E+C1       0.1894E+03  0.6563E+01  0.3020E+00 

25  12            21       0.5979E+C2       C.10S1E+03  C.19C0E+02  0.1002E+01 

26  11             IS       0.5979E+02       0.1081E+03  C.1900E+02  0.10C2E+01 

27  9            16       C.5979E+C2       0.1111E+03  C.1900E+02  C.1302E+01 

28  10            17       C.5979E+C2       0.1111E+03  C.  lQ00E+02  G.1CC2E+01 

29  9            21       0.7358E+C1       0.2380E+03  0.77o7E+01  G.3C16E+00 
3C          1C            16       C.73?eE+Cl       0.2380E+03  C.7767E+01  0.3016E+00 

31  11            13       0.7358E+Q1       0.6846E+02  C.7767E+01  0.3C165+00 

32  12            13      0.7358E+01       0.6846E+02  C.7767E+01  0.3016E+00 
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33 

13 

18 

G.7256E+C1 

34 

13 

21 

C.73  5  6E+C1 

35 

1C 

Its 

0.7358E+01 

36 

16 

21 

C.  5424E  +  C1 

37 

17 

18 

0.5424E+C1 

38 

18 

21 

G.5424E+01 

3<5 

16 

17 

0. 5424E+C1 

*0 

12 

20 

0.5424E+G1 

41 

11 

IS 

C.5424E+C1 

42 

20 

21 

C.5424E+C1 

43 

18 

19 

0.5424E+01 

44 

19 

2C 

0. 5424E+01 

45 

21 

29 

0.597<;5  +  C2 

46 

18 

28 

O.5979E+02 

47 

16 

26 

0.5979E+C2 

48 

17 

27 

C.5979E+C2 

49 

16 

29 

C.10C8E+C2 

50 

17 

26 

0. 1C<SE+C2 

51 

20 

22 

Q.5424E+01 

c  5 

19 

2C 

C.5424E+C1 

53 

20 

21 

O.1098E+G2 

54 

'     19 

al 

0.10985+02 

55 

18 

24 

0.10*85+02 

56 

21 

24 

C.1G9SE  +  C2 

57 

24 

26 

C.  1QS8E+C2 

58 

24 

29 

0.  1098E+02 

59 

16 

25 

C.  1C<8E+G2 

60 

17 

25 

0.1098E+O2 

61 

26 

29 

0.  1648E+02 

62 

27 

2S 

0.  1646E+C2 

63 

29 

22 

0.  1648E+02 

64 

28 

2C 

C.  16A8E  +  32 

65 

21 

22 

0.  1648E  +  C2 

66 

30 

21 

0.1648E+02 

67 

29 

29 

0.1646E+C2 

68 

15 

6C 

0.2142E-C6 

69 

15 

6C 

0.2142E-06 

7C 

2C 

22 

C.  16^6E  +  C2 

71 

29 

29 

0. 5979E+C2 

72 

28 

27 

0.5979E  +  C2 

73 

26 

24 

C.59795+C2 

74 

27 

26 

0.5979E+C2 

75 

26 

35 

0.1667E+C2 

76 

27 

27 

C. 1667E+C2 

77 

29 

40 

0.1667E+02 

78 

2e 

26 

C.  16675  +  02 

79 

26 

52 

0. 1667E+02 

80 

27 

52 

0.  1667E  +  C2 

81 

28 

~.  -a 

0.  1667E  +  C2 

82 

29 

33 

0.  16675+02 

93 

33 

27 

0. Ic675  +  C2 

84 

■a -a 

-c 

0.1667E+C2 

85 

34 

39 

0. 16675+02 

ge 

26 

27 

C.  16675+C2 

37 

39 

*G 

0. 1667E+G2 

38 

37 

23 

0.  1667E+C2 

39 

39 

47 

C.  5979E  +  C2 

90 

37 

45 

0.5979E+02 

91 

24 

t-2 

0.597CE  +  C2 

92 

36 

i? 

C. 5C79E+C2 

93 

34 

47 

0. 16975+02 

94 

26 

*5 

C.  1697E+C2 

95 

40 

*8 

0.  1667E  +  02 

96 

38 

46 

0.  1667E+02 

97 

25 

It 

G.  1667E+C2 

98 

25 

27 

0.16675+02 

99 

26 

1  c 

C.16675+C2 

100 

27 

-.  c 

0.  16675+02 

101 

24 

35 

0.1667E+02 

102 

■3  C 

26 

C.  1667E  +  02 

10  3 

35 

42 

0.  1667E  +  02 

1C4 

35 

43 

0.1667E+02 

0.68465+02  0.7767E+01  0.3016E+00 

0.6846E  +  02  C.7767E+01  0.3C165+00 

0.2094E+03  0.7767E+01  0. 3016E+OO 

0.18665+03  0.50245+01  0.3437E+0O 

C. 18666+03  C.  50245+01  C. 34375+00 

0.34C0E+02  0.5024E+01  0.3437E+00 

C.1663E+03  C. 50245+01  0.3437E+00 

0.1222E+03  C.5024E+01  0.2427E+OO 

0.1232E+03  G.5024E+01  0.3437E+00 

0.6456E+02  C. 50246+01  0.34375+00 

0.64566+02  0.5O24E+01  0.3437E+00 

0.84COE+02  0.50245+01  0.3437E+O0 

0.1081E+03  C.19CCE+02  C.10C2E+01 

0.1031E+03  C.19C0E+02  0.10026+01 

0.1112E+03  C.1900E+02  0.1002E+01 

0.11125+0?  C.19CCE+02  G.1CC2F+01 

0.2112E+03  C.9952E+01  0.35115+00 

C. 21125+03  C.9952E+01  0.3511E+00 

0. 12095+03  C.5024E+01  0.3427E+00 

C.12095+03  0.5024E+01  0.3437E+00 

0.1626E+03  C.9952E+C1  0.2511E+O0 

0.1626E+03  C. 99525+01  0.35115+00 

0.68*65+02  C.9952E+C1  0.2511E+OO 

0.68465+02  0.9952E+01  O.3511E+0O 

C.6846E+02  C.  9952E+01  0.2511E+O0 

0.6846E+02  C.9952E+01  0.25UE  +  00 

0.1334E+03  G. 99525+01  0.3511E+00 

0.13846+03  C.9952E+01  C.3511E+00 

0.1558E+03  0.10825+02  0.48485+00 

G.  1558E+03  C.i0325+02  0.4848E+00 

Q.1240E+02  C. 10825  +  02  0.4848E+00 

0.1240E+03  0.10825*02  0.4348E+00 

0.76845  +  02  0.1082E  +  02  0.4848E+O0 

0.768*6+02  C.1032E+02  0.4848E+00 

0.840CE+02  0.1082E+02  0.48485+00 

0.10COE+00  CO  0.0 

0.1000E+00  0.0  0.0 

C.96C0E+02  C.1082E+02  0.484SE+00 

0.96115+02  C.19C0E  +  02  C.1GC2E+01 

0.9611E+02  0.19005+02  C.10G25+01 

0.9979E+02  C.1900E+02  0.10025+01 

0.9379E+02  C.1900E+02  C.10C2E+01 

G.1792E+03  0.1056E+02  0.5O26E+OO 

0.1792E+03  C.1056E+02  0.5026E+00 

0.10005+03  0.1056E+02  0.50265+00 

G.  10005+03  C.1056S+C2  0.5026E+00 

0.20665+03  C. 10565  +  02  0. 30265+00 

C. 20665+03  0.10565+02  0.5026E+00 

0.63825+02  C.1056E+02  0.5026E+00 

0.6332E+02  O.I056E+02  0.5O26E+0O 

0.6382E+02  0.10565+02  0.50265+00 

0.6382E+02  C.  10565+02  0.50265+00 

0.12895+03  0.10565+02  0.5026E+00 

C.12S9E+03  C.1056E+02  0.50266+00 

0.32545  +  C2  C.  10565+02  0.5026E+OO 

0.32545+02  0. 10565+02  0.50265+00 

0.  72C9E+02  C.1900E+02  0.10025+01 

0.7209E+02  0.1900E+02  0.10C2E+01 

0.74116+02  0.1900E+02  0.10025+01 

0.74115+02  C.1900E+02  C.10C2E+01 

0.1446E+03  0. 10355+02  0.5219E+OO 

0.14465  +  03  C.1035E+02  0.52195+00 

0.7428E+02  C.1056E+02  0.5C26E+00 

0.74285+02  0.1056E+02  0.5026E+00 

C.7176E+02  C.1056E  +  02  0.5026E+00 

0.7176E+02  C.1056E+02  0.5026E+00 

0.1217E+O3  0.10565+02  0.50266+00 

0.12175+03  C.  10565+02  0.5C265+00 

0.61685+02  0.10565+02  0.5G26E+00 

0.61&65+02  C. 10565+02  0.5026E+00 

Q.9138E+02  C.1056E  +  02  O.5026E+0O 

C.9138E+02  0.10565+02  0.5026E+00 
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VII .   CONCLUSIONS  AND  RECOMMENDATIONS 

A.   CONCLUSIONS 

An  existing  finite  element  code  was  expanded  to  encom- 
pass the  more  general  case  of  frame  members;  i.e.,  six 
degrees  of  freedom  per  joint.   Combined  truss  and  frame 
structures  were  designed  for  minimum  weight  with  multiple 
load  conditions  considered. 

The  displacement  method  was  used  for  static  analysis 
and  the  subspace  iteration  method  was  used  for  eigenvalues. 

Several  examples  were  considered.   In  every  case  the 
code  worked  as  an  analysis  tool,  and  significant  weight 
reductions  were  obtained  with  the  coupled  optimizer  CONMIN. 

The  SADX  design  code  has  been  shown  to  be  a  useful  tool 
for  ship  mast  optimum  design. 

3 .   RECOMMENDATIONS 

The  following  recommendations  may  be  of  value  for  future 
work . 

1.  The  routines  necessary  to  calculate  gradients  analyt- 
ically  should  be  added  to  the  code. 

2 .  The  code  should  be  extended  to  include  plate  and 
shear  elements  and  a  library  of  member  cross-sections. 

3.  An  out  of  core  ecuation  solver  should  be  added. 
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M- .   _he  ability  tc  specify  multipliers  for  applying  ir.er- 
tial  loads  along  the  three  coordinate  axes  should  be  added 
in  a  fashion  similar  to  that  used  for  applying  structure's 
own  weight  as  loads.   Such  an  addition  would  simplify  dynamic 
load  analysis  and  design. 

5.  The  method  of  gradient  calculation  should  be  depen- 
dent en  specific  gradients  recuired  [Ref .  7  ]  and  [Re:  .  z~\  . 

6 .  Gradients  cf  frequency  const: 
a  tore  efficient  ai  ^or i —  ht.  LRef .  3  J  . 

7.  The  need  for  a  larse  scale  tub  lie  structural  cttir.iza- 
ticn  code  still  exists . 

3.   The  code  should  he  modified  to  allow  optimum  member 
size  desizn  followed  by  simultaneous  Dptimum  member  size  end 
optimum  geometry  design. 
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APPENDIX  A 
DATA  FILES 

A.   INTRODUCTION 

This  appendix  contains  the  data  files  used  to  create 
the  test  cases  in  Chapter  VI.   Additionally  rhe  data  file 
for  the  USER'S  guide  in  complete  form  is  presented. 
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TABLE  X 
DATA  FILE  TRUSS-BRACED  CANTILEVER  BEAM 

$*«*** 4*41***4 4****2  ******** *3 ******  ***^** ******* 5** *******£ ********* 7********* 
SBLCCK     A    TITLE  1    CARD  ^ORMAT       30H 

1RLSS-BRACSD    CANTILEVER    BEAM:    EXAMPLE    1 

$*****  *  *  *!*:**  44*  4**2  4*4******3*********^**  *******  5**  **  *****£***  ******  7**  ****  *** 

SBLCCK  B  CCNTPCL  PARAMETERS 

J   NEB        NEF         NJ        NCJ 

2  2  5 

S  NEUBC      L8LCK   ,   NFREQ   .  NFMASS 


I 


3 

1 

CARDS 
NMT 

FORMAT 
IOVCLC 

3110 
■ 

NDJ 

I 

i 
3 

1 
1 

l 

NEIG 
NFMW 

2 

NEIG1 

l 

2 

1 

NPRI 

2 

i 

i 

1 

! 
1 

$      NLC  NCSPLC  NSTR5S  NSTW 

Sill 
2  2  2 

J* 4*** 4  44 1*** 4 < 4***2 4**4***** 3  ****** 4**4**  4*** *** 5********* A**** *****  7** **** *** 

S9LGCK    C    DYNAMIC    ANALYSIS     INFORMATION       1    CARD  FORMAT       I10,2Fi0.0 


.1  -' 


0  386.4  0. 

«**4***4 4  1** 444****24*44*****3*********4** *******5*********6*4*******7********* 

SBLCCK    C    JOINT    COORDINATES 


$            NJ    CAPOS 

FOPMAT 

I10.3F10 

$         JN              X-CCCRO         Y-COORO         Z-CCORO 
c                        I                           1                           1                           I 

*                                               1                         1 

0.                       0.                       0. 

2                  100.                       0.                       0. 

3                 200.                      0.                      0. 

4                      0.                  150.                       0. 

5                      0.                       0.                     50. 

SBLCCK     E    DESIGN    VARIABLE    LINKING    DATA 

i            NDJ    CAROS 

FGR"*"r 

4I10t3Fl0.0 

$          JN                     IX                        IY                        IZ 

PCX 

PCY 

PCZ 

Still 

4                         C                         1                         0 

I, 

1 
.0 

1               1 

1.0                    1.0 

5                         0                         0                         2 

1. 

.0 

1.0                    1.0 

$********! ********  * 2 ******** *3 ******  ***4**  **** ***5*******:**64********7*** ****** 
SBLOCK    F    MftTERIAL    PROPERTIES    "  NMT    CAPOS         FORMAT       6F10.0 

$  E  RhC  SIGMIN  SIGMAX         KEULER  POISSN 

1.0E+7                   0.1  -25000.            25C00.  2.                    .27 

2.SE+7                    0.3  -36000.             36000.  2.                      .27 
J********l***44****2 4** 4*****3  ********* 4**  4***  *** 5** ******* 6***** 4*** 7* ******** 

SBLCCK    G    EAR     ELEMENT  INFORMATION               NEB    CARDS  FORMAT       5I10fF10.H0 

S         LNC              NCCE2  NCDE3            MATCOD               NSDG1               AREA                 LSECT 

5                        I                            I  I                            I  I                            I                           I 

1  2  4                         1  1                    3.0                         1 

2  2  5                         1  2                    3.0                         1 

$*****4441***»4****244**4* 4* 43 4********4**  4****** 5 *»*******£********* 7********* 

SBLCCK    h    FPAVE    ELEMENT     INFORMATION  2  *N  EF    CARDS     FORMAT    7I10/2F10 

S         LNO  N0CE2  NCDE3  MATCOD  NS0G1  NSDG2  LSECT 

5  I  I  I  I  I  I  I 

3  12  2  3  4  1 
ICI-AR0IV1      CHAPCIM2 

S  I  I 

;.c  1.0 

4  2  3  2  3  5  1 

5.C  1.0 

******  44*1***44***  *2**4******'a*********^**  *******5  S********A ********* 7 ********* 

SBLOCK     I    JOIN'    CONSTRAINT    DATA  NCJ    CARDS  PORMAT       7110 

S       JN  IX  IY  IZ  IXX  IYY  IZZ 

I  1  1  i  1  I  1 

2  0  0  0  0  0  0 

3  0  0  0  0  0  0 

4  111111 

5  1111  11 

J* 4*****4 1*4* ******  24  44  4*****3  4***4****4**  4*** *** 5** ** **** *6*** * *****  7****** *** 

SBLCCK    J    JCIMT    LGACI NG    CATA  1    CAPO  FORMAT     110 

i  <Cf!7    IF    NLC=C    IN    8LCCK    B) 

S       NLJ 

S  I 

1 


77 


TABLE  XI 
DATA  FILE  TRUSS-3RACED  CANTILEVER  3EAM  continued 

*  ,„  _„  m  NLJ    CAROS  FORMATI10,6F10 

i     '"     i       FX     i       "     .       FZ     |       TX     i       "     i       Ti     | 

3  o,         laco.o  o.  o.  o.  o. 

3  0.  0.0         -1000.0  0.  0.  0. 

j*  *******!********;:*  2*********3  ******  ***4»*  **** 4**5**** 9** **£******* 44 7********* 

SBLCCK    K    FIXED    MASS     CATA  NFMASS    CARDS         cORMAT    U0,F10.0 

$  (CKIT    IF    NFMASS  =  0    IN    BLOCK    B) 

J         JN  MASS 

$  I  I 

3  250. C 

j*  ******  *1***  *****  *?*********  3*  *******4<u**  44*****5 *********£ ********* 7 ********* 
SBLCCK    I    DESIGN    VARIAPLF     INFORMATION       CARDS    AS    REQ'D  FORMAT    3F10.0 

S  (AREA    ANC    CIMENSICN    VARIABLES  CMIT    IF    NDVAP1  =  0) 

$  XA(I)    .        .       .        XA(NOVARl) 

$11  I                      I                   I                    I 

2C.C                20. 0  20.0                   2.0                   2.0 
S         XAL(I) 

$                       I                          I  I                          I                            III 

.50                  .50  4.0                   .10                     .10 
S         XAU(I) 

i  '   I  1  I  J  ill 

25.0  35. C  25.0  2.5  4.0 

<:s4****  4*1*********2  ********* 3 ********* 4**  ******* 5** ******* 6 ********* 7********* 
SBLCCK    ¥    DESIGN    VARIABLE     INFORMATION       CARDS    AS     REO'D  FORMAT    8F10.0 

J  [COORDINATE    VARIABLES  OMIT    IF    NDVAR2=0) 

S  XCm    .        .      .       XC(NDVAR2) 

sill  III 

150.0  50.0 

S         XCL(I) 

s  I  I  I  I  III 

0.0  0.0 

sill  111 

200.0  1C0.0 

«********!**»*  ***** 2*********3  *********<►**  *******e*********6 ********* 7********* 

SBLCCK     !\    JCINT    C!S?LJCCMENT    CONSTRAINTS  NOSPLC    CAROS       FORMAT    3U0t2F10 

S  (CMIT     IF    NDSPLC=C     IN    BLOCK    3) 

i         JN  CIR  LC  SL  BU 

$  I  I  I  I 

3  2  1-3.0  3.0 

3  3  2-3.5  3.5 

5********1  *********•?****** ***3  444*4* ***A**4******5  ********* 6*********7********* 
iBLCCK    C    FREQUENCY    CONSTRAINTS  NO.    OF    CARDS    AS    REQ'O  FCRMAT    3F10.0 

«  (CVIT     IF    NREC=0     IN    BLOCK    3) 

S  FC1       .  FCN 

still  III 

1.0 
j******** i** ******* 2*********3  ******  ***4**  *******5*********6********* 7 ********* 
SELCCK    P    ENO    C*RC 
ENC 
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TABLE  XII 
DATA  FILE  TWO-TIER  3-D  PORTAL  FRAME 


I******** l** ****** 42*********3  ********* 4**  *******5*********6 ********* 7 ********* 

jeiCCK    A    TITLE 

1    CARD            FORMAT       80H 

2-TIEF    FRAM  =    STRUCTURE    WITH    TRUSS    X-BRACcS:    EXAMPLE    2 

$4* ******  I** ******  a 2*********3 ******  ***4**  ******* 5** »4*****6*********7 ********* 

JEICCK     E    CONTRCL    PARAMETERS                                  3    CARDS         FORMAT       9110 

i       NEB                     N-F 

NJ                    NCJ                    NMT               IDVCLC                 NO J 

i                 1                   1 

20                       16 

il             J             1             i             1 

$    NEUBC               L8LCK 

NFREQ            NFMASS              NE IG                 NEIG1              NPRI 

1             1 

1             i             1             I             I 

f       NLC              NOSFLC 

NSTRSS            NSTW                 NFMW 

3                       12 

\           !           1 

J********]*4  4  4*****2*********'3  4***4****4**  4******5*******4*6*********7********* 

SELCCK    C    CYNAMIC    ANALYSIS     INFORMATION       1    CARD            FORMAT       I1C2F10.0 

$    LMASS              GR/V 
$                       1                          1 
0                3  S6  .  4. 

EPSEIG 

l 

0. 

«* *44* ** *i ** ****** 42  ***** ****3 ****** ***4*********5**** ***** 6*** ******  7****4**** 

J3L0CK    C    JOINT    COORDINATES 

$            NJ    CARDS 

FORMAT    I10,3F10 

$         JN              X-CCORO 

1                      i                          > 
1                 100. 

Y-COORD         Z-COORO 

i                          i 

1                          1 
0.                 100. 

2                 100. 

0.               -100. 

3              -100. 

0.               -100. 

4              -100. 

0.                  100. 

5                 100. 

50.                  100. 

6                 10  0. 

50.               -100. 

7              -100. 

50.               -100. 

8              -IOC. 

50.                  100. 

9                 100. 

100.                  100. 

10                  100. 

100.               -100. 

11              -100. 

100.               -100. 

12               -100. 

100.                  100. 

$***** 4** 1*4  4 ****** 2 ******** 41  ******  4**4**  ******* 5**** 4****6********* 7********* 

JELCCK     E    CESIGN    VARIABLE    LINKING    DATA 

S            NCJ    CARDS 

FORMAT    A 1 10, 3F  10.0 

1          JN                     I) 

{              1 

IY                      IZ                    PCX                   PCY                    PCZ 

0                         2                    1.0                    0.0                    1.0 

2                        1 

0                         2                    1.0                    0.0                 -1.0 

3                        1 

0                         2                 -1.0                    CO                 -1.0 

4                         1 

0                         2                 -1.0                    0.0                    1.0 

5                        2 

4                         5                    1.0                     1.0                    1.0 

6                        3 

4                         5                    1.0                     1.0                 -1.0 

7                        3 

4                         5                 -1.0                     1.0                 -1.0 

8                         3 

4                         5                 -1.0                    1.0                    1.0 

«.*****  ***]  *********2  *********?  *********  A**  *******  5  *4****4  4  4  6*****44  4*  7********* 

1BL0CK    F    MATERIAL    PRCPERTIES    "                      NMT    CARDS          FORMAT       6F10.0 

S            E                    RJ-C 

$               I                 ! 

2.9E+7                   0.3 

SIGMIN             SIGMAX         KEULER            POISSN 
-A2C00.             36C00.                        4.                      .27 

J*  44****4 1*******4*2** 4*4*4**3  ********* A** ******* 5*********6*** ******  7********* 

53LCCK    G    BAR     ELEMENT 

INFORMATION               NEB    CARDS          FORMAT       5110, F10, 110 

i         LNC              N0CE2 

NCDE3            MATCOC               NSDG1               AREA                  LSECT 

1             1 

I                       1                       t                  5.i 

2,5,2,1,1,5.0,1 

3,2,7,1  ,1,5.0,1 

A, 6, 3, 1,1, 5. 0,  1 

5, 3,9,1,1,5.0,  1 

6,7,4,1  ,1,5.0,1 

7, A, 5,  1,  1,5.0,  1 

6,1,8,1,1,5.0, 1 

C5, 7, 1,1, 5. 0,1 

10,6,8,  1,1,5.0,1 

11, 5, 1C, 1,2, 5. CI 

12, 6,  =  ,1,2,  5.  CI 

13  ,6, 11, 1,2, 5. C,i 

14, 7, 1C, 1,2,5. CI 

79 


TABLE  XIII 
DATA  FILE  TWO-TIER  3-D  PORTAL  FRAME  continued 

15, 7.12,1,2,5.01 

16, 8, 11, 1,2, 5. C.l 

17,8,9,1,2,5.0  ,1 

IP, 5, 12, 1,2, 5. C,l 

19, 9, 11, 1,2, 5. C,l 

20,10,12,1,2,5.0,1 

i  ***=***»*  I**  *******  2*********3*********4**  *******  5**  *******6*********7********* 

S3LCCK    N    FRAME    ELEMENT     INFORMATION  2*NEF    CARDS    <=ORMAT    7I10/2F10 

S         LNG  NCCE2  NCDE3  MATCOC  NSDG1  NS0G2  LSECT 

Zl  tl  ,5*1,3,4*1 

5.C  0.50 

22,2,6,1,3,4,  1 

5.0  0.50 

22,3,7,1,3,4,1 

5.0  0.50 

24  ,4,8,1,3,4,1 

5.C  0.50 

25,5,9,1,5,6,1 

5.0  0.50 

26,6,10,1,5,6,1 

5.0  0.50 

27.7.11  ,1,5,6,1 

5.C  C.50 

28.8.12  ,1,5, 6, 1 

5.0  *  0.50 

29,5,6,1,7,8,1 

5.0  0.50 

30,6,7,1,7,8,1 

5.C  C5G 

31,7,3,1,7,8,1 

5.0  0.50 

22,6,5,1,7,8,1 

5.0  0.50 

22, 9, 1C, 1,9,10,1 

5.0  C.50 

34, 10, 11,1, 9, 1C,1 

5.C  0.50 

35,11, 12,1,9, 1C,1 

5.0  0.50 

26, 12, <=, 1,9,10.1 

5.C  C.  5C 

<***»*:»**]_*«***»***?***  ******3  ********* 4**  4  ***  ***  5  *********  £********* 7***  ****** 

J3  1CCX     I    JOINT    CCNSTPAINT    DATA  NCJ    CARDS  FORMAT       7110 

J       JN  IX  IY  IZ  I XX  IYY  IZZ 

J  I  I  I  I  I  I  I 

1,1.1,1,1,1,1 

2,  1,1,  1,1,1,1 

3,1,1,1,1,1,1 

A, 1,1, 1,1, 1,1 

5, 

6, 

I: 

9, 

10, 

11, 

12  , 

j*****4**l***«  <**** 2*********3*********^** *******5*********6 *********  7 ********* 
JBICCK    J    JOINT    LCACING    CATfl  1    CARD  FORMAT    110 

i  ICflT    IF    M_C=0     IN    3LGCK    B) 

1      NLJ 

I 

NLJ   CARDS 
FX  FY  FZ  TX 


J 

1 

2, 

S 

i 

JN 

$ 

1 

9, 

1000. 

1C 

,1000. 

«, 

«, 

500. 

10 

,500. 

FORMAT 

U0,6F10 

TY 

TZ 

1 

1 
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TABLE  XIV 

DATA  FILE  TWO-TIER  3-D  PORTAL  FRAME  continued 

5t-5CO. 
t-5CC. 

C,1000. 

.C.CC.C1000. 

2,-lOCC. 

1  »C.C,O.C,-10CO. 

********  i**  *******  2*********3  ******  ***4.**  4»«*««5»*$4«4«^»«****«**7«******** 
S3LCCK    <    FIXEC    "ASS    CATA  NFMASS    CARDS  FORMAT     I1CF10.0 

5  <CMT     IF     NFMASS=0     IN    3LOCK    3) 

JN  M A <  < 

,50.0 

0.50.C 
1,50.0 

2  ,5C.C 

******  **!***  ******  2*********3  ******  *******  ******* 5** ******* a ********* 7** ******* 

eLCCK    L    CESIGN    VARIA31E    INFORMATION      CARDS    AS    REQ'Q  =CRMAT    8F10.0 

(AREA    ANC    CIMENSICN    VARIABLES  CMIT     IF    NCVAR1  =  0) 

XA(I|    ...        XAtNOVARll 

5  3.5  5.0  0.5  5.0  0.5  5.0 


XAK  II 


2.5  3.5  5.0  0.5  5.0  0.5  5.6  C.5 

e.O  0.5 


1.1  .50  0.1  .56 


0.1                   0.1                   .50  0.1  .56                 0.1                   .50                0.1 
.50                   0.1 
XAL(I) 

I                          I                          I  I  III 

3C.0                 30.0                 30.0  .50  30.0                  .50                 30. C                 C.5 
2C.C                   . 5C 

******4*£****  4*  ***2*********'a  ******  ****.**  *****************  a***  ******  7******  *** 

eLCCK    F    DESIGN    VARIABLE     INFORMATION  CAROS  AS     REQ'D             FORMAT    3F10.0 

(CCCRDINATE     VARIAELES  CMIT  IF     NDVAR2=0J 


XC(I)  .        .       .        XC(N0VAR2) 
I  I  I 

1CC.0  1C0.C  1C0.O  ICO. 0 

XCKI 


>o.i 


I 


III  III 

2C.0                 20.0                 20. 0                 20.0  20.0 

xcum 

Fill  ill 

16C.0                160.0                160.0               160.0  30.0 
******  *  *  i  ********* 2  ***  ******?  ******  ***4**  *******=  *********«.  ********* -7  ******  *** 

6LCCK    N    JOINT    DISPLACEMENT    CONSTRAINTS  NDSPLC    CARDS       FORMAT    3110, 2F10 
(CUT     IF    NCSPLOO    IN    3L0CK    3) 

JN                    OIR                    LC                         8L  3U 

-5.0  5.0 

-5.0  5.0 

-5.0  5.0 

-5.0  5.0 

-5.0  5.0 

-5.0  5.0 

-5.0  5.0 

-5.0  5.0 

-5.0  5.0 

-5.0  5.0 

-5.0  5.0 

-5.0  5.0 
j********  1****4**** 2 ********* 3*=****  *** 4**  *******5 ****=*****,********* 7 ********* 

SELCCK    C    FREQUENCY     CONSTRAINTS             NO.    OF  CARDS    AS    REQ'D             cORMAT    3F10.0 

S  (CVTT    IF    NFREC=0     IN    SLCCK    3) 

J  FC1       .       .    FCN 

S                       I                           I                           I                           I  III 

l.C  5.0  10.0 

4* :»  *****<£*********  2*********3******  ***4**»******5*********6*********  7********* 
5ELCCK    F    END    CXRC 
ENC 


9 

1 

1 

10 

1 

1 

11 

1 

1 

12 

1 

1 

S 

1 

2 

10 

1 

2 

11 

1 

2 

12 

1 

2 

9 

3 

3 

10 

1 

3 

11 

3 

3 

12 

1 

3 

Bl 


TABLE    XV 
DATA   FILE    DD-96  3    FOREMAST 

SBLOCK    A    TITLE 

CC-963    CLASS    CESTROYER    FORWARD    MAST    REDESIGN:    EXAMPLE    1 

$       THIS     IS    A     SINPLIPICiTICN       GF    THE    ACTUAL     STRUCTURE    ANO    LCAOS    8ASEO    QN 


S       NAVSHIPS 

DPAWING    MiM 

EER.    128- 

453S510    AND    NAVSHIPS    SKETCH    NUMBER 

J       8006-1-123 

-SK45379  79 

$       THIS    ANALYSIS    RECLIPES     A    THE 

X-AXIS    FORWARD,    Y-AXIS    UP,    AND    Z-AX  IS 

*       TO    ST60; 

A    CIFFERENT 

CRIENTATION    THAN    IN    THE    ABOVE    DRAWINGS 

SEICCK    6    CCNTPCL    PARAMETERS 

i       NEB 

i                       1 
0 

NEF 

1 

NJ 

NCJ                    NMT               IDVCLC                 NDJ 

1                          1                          1                          I 
20                        1                        1                        0 

172 

73 

t    NEUBC 
S                        1 
1 

LBCCK              NFREQ 

NFMASS              NSIG                 NEIG1              NPRI 

1 

1 

0 

33                         0                         0                         0 

$       NLC              NDSFLC              N 
<                       I                          * 

STRES 

■ 

NSTW                 NFMW 

i                          i 

1 

ll 

2 

1                         1 

$*** *****!** 

*******  2  **=****:***3  ******  ***^**  ******* 5*********6********* 7*** 

SELCCK    C    CYNAMIC    ANALYS 

IS     INFORMATION 

C 

386.4 

0. 

JBLOCK    C    JOINT    CCOROINAIE 

1 

-240.0 

0.0 

-1C6.0 

2 

-240.0 

0.0 

1C6.0 

a 

O.C 

0.0 

42.0 

4 

0.0 

0.0 

-42.0 

5 

-223. C 

54.0 

0.0 

6 

-120. C 

54.0 

6  8.4 

7 

-2.64 

54.0 

0.0 

8 

-120. C 

54.0 

-68.4 

9 

-216.C 

1C8.0 

-«4.  7 

10 

-216.0 

1C3.0 

94.7 

11 

-5.  28 

1C8.0 

42.0 

12 

-5.28 

1C8.0 

-42.0 

13 

-7.92 

162.0 

0.0 

14 

-172. C 

444.0 

0.0 

S    NOCE    13    IS 

A^    UNCCNNECTEC    FIXED    DUMMY   NODE 

15 

0.1 

0.0 

0.0 

16 

-192.6 

216.0 

-83.16 

17 

-192.6 

216.0 

83.16 

18 

-1C.56 

216.0 

42.0 

19 

54.0 

216.0 

42.0 

20 

54.0 

216.0 

-42.0 

21 

-1C.  56 

216.0 

-42.  0 

22 

-172  .C 

726.0 

-18.0 

22 

-172.0 

736.0 

18.0 

24 

-13.2 

27G.0 

0.  0 

25 

-168.8 

324.0 

0.0 

26 

-163. S 

324.0 

-71.76 

27 

-168.8 

224.0 

71.76 

28 

-15.84 

324.0 

42.0 

29 

-15.  84 

224.0 

-42.0 

20 

1C8.0 

224.0 

48.0 

21 

168  .0 

224.0 

0.0 

22 

1C8.C 

224.0 

-48.0 

23 

-18.  12 

272.0 

0.0 

34 

-147 .8 

420.0 

-6  1.6  8 

:e 

-147.8 

420.0 

0.0 

26 

-147.8 

420.0 

61.63 

27 

-20.4 

420.0 

42.0 

28 

12. C 

420.0 

29.0 

29 

-20.4 

420  .0 

-42.0 

4C 

12.0 

420.0 

-29.0 

41 

-22.2 

456.0 

0.0 

42 

-122.0 

492.0 

-54.0 

43 

-132.  C 

492.0 

54.0 

*4 

-172.0 

492.0 

-54.0 

45 

-24. C 

492.0 

42.0 

46 

20.0 

492.0 

36.0 

47 

-24.0 

492.0 

-42.0 

4E 

20.  C 

492.0 

-36.0 

$    ENC    NCCES 

FCF    THE    LOWER    QUADRAPOD    CF    THE     FORWARD    MAST 

S    START    THE 

UPPER    FCFEfAST    NOOES 
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TABLE  XVI 
DATA  FILE  BD-9  6  3  FOREMAST  continued 


i    NODE 


$****** 

1SLQCK 
2.9 

iBLOCK 
5  INC 

i 
******* 

setccK 

5  START 
11 
11 
11 
11 
7. 
7. 
7. 
7. 
7. 
7. 
7. 
7. 
7. 
7. 
7. 
7. 
7. 
7. 


49 

50 
51 
52 

53 

54 

ec 

56 
57 
58 
59 
6C    IS 

to 

61 
62 
63 
64 
65 
66 
67 
68 
69 
"iC 
71 
72 
72 
**!**** 

F    PATER 
E  +  7 
**!**.** 

G    eAR    E 
NC 
I 
**]**** 

h    FRAME 
FRAME 

1 
.75 

2 
.75 

•a 

.75 

4 

.75 

c 

625 

6 
625 

7 
625 

e 

625 

9 

625 

10 
625 

11 
625 

12 
625 

13 
625 

14 
625 

15 
625 

16 
625 

17 
625 

13 
625 


-172. C 

-172. C 
-228.0 
-2?3.C 
-2S8.C 
-2^0.0 
-172. C 
-172.0 
-213.0 
-213.0 
-172.0 
N    L'NCCNNE 

0.0 
-172.0 
-172.0 

-  1  72  .  C 
-172.0 
-172.  C 
-172.0 
-172.0 

-  1  72  .  C 
-172.0 
-172. C 
-172.0 
-172.0 
-172.C 

t  >  »<:»2  *a  * 
1AL     INFCP 

0.1 
*****2*** 
LEMENT    IN 
CE2  N 

I 
*****2**< 

ELEMENT 
ElEMEMS 

4 
.75 

3 
.75 

1 
.75 

2 
.75 

4 
.375 

3 
.375 

1 
.3  75 

2 
.375 

8 
.375 

6 
.375 

8 
.3  75 

6 
.375 

3 
.375 

4 
.375 

7 
.375 

7 
.3  75 

1 
.3  75 

2 
.375 


492.0  -120.0 

492.0  120.0 

408.0  0.0 

492.0  -56.0 

492.0  56.0 

492.0  0.0 

492.0  -276.0 

492.0  276.0 

600.0  -24.0 

6C0.0  24.0 

492.0  54.0 
CTEC    FIXED    CUNMY    NOCE 

0.0  0.0 

600.0  0.0 

492.0  0.0 

664.0  0.0 

736.0  0.0 

396.0  0.0 

976.0  0.0 

1056.0  0.0 

736.0  -264.0 

736.0  264.0 

736.0  -156.0 

736.0  156.0 

736.0  -228.0 

736.0  228.0 
»***433*****s***4***«*****5  *********£* =*******  7*** 

f-ATICN 

-42000.  42C00.  4.  .27 

*****  *3  *********<,**  ******* 5 *********£********* 7*** 

FORMATION  NEB    CARCS  FORMAT       5110,  F10.U0 

CDE3  MATCOC  NSCG1  AREA  LSECT 

I  I  I  I  I 

****** 3 ********* 4* ********5  ******  ***£,********* 7*** 

INFORMATION 

FCR    1ST    CR    LOWEST    TIER    OF    QUAQRAPOC 

1  2 


*3 
MA 
ST 
12 

11 
9 

10 
3 

6 
8 
6 

12 

11 

9 

10 
7 
7 
11 
12 
5 
5 
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'ABLE   XVI 


DATA   FILE    22-3H    FOREMAST    continued 


19 

5 

<?                          1 

7.625 

.3  75 

2C 

e 

10                          1 

7.425 

.375 

21 

12 

9                         1 

7.425 

.375 

22 

11 

10                        1 

7. £25 

.375 

23 

11 

12                        1 

7.625 

.375 

24 

9 

10                       1 

7.625 

.375 

J    EMC    1ST   T!=b 

CP    CUACR4PQC    START    2N  C    TI 

25 

12 

21                        1 

11.75 

.75 

26 

11 

ia                 l 

11.75 

.  75 

27 

9 

16                        1 

11.75 

.75 

28 

1C 

17                          1 

11.75 

.75 

2  = 

9 

21                         1 

7.625 

.375 

30 

10 

18                          1 

7.625 

.375 

31 

11 

13                        1 

7.625 

.3  75 

32 

12 

13                        1 

7.625 

.375 

■=■5 

13 

13                         1 

7.625 

.375 

34 

13 

21                        1 

7.625 

.3  75 

■:  e 

1C 

16                         1 

7.625 

.375 

36 

21 

16                         1 

7.625 

.375 

37 

17 

13                         1 

7.625 

.3  75 

3  8 

18 

21                         1 

7.625 

.375 

3« 

16 

17                         1 

7.625 

.375 

40 

12 

20                         1 

11.625 

.3  75 

-1 

11 

19                         1 

11.625 

.375 

*2 

20 

21                         1 

7.625 

.375 

4"= 

13 

19                          1 

7.625 

.375 

44 

19 

20                         1 

7.62  = 

.375 

5    ENC    THE    2N0 

TIER     CF    THE     CUAORAPGC    STAR- 

*5 

21 

29                         1 

11.75 

.75 

46 

IS 

23                        1 

11.75 

.75 

47 

16 

26                        1 

11.75 

.75 

43 

17 

27                       1 

11.75 

.  75 

49 

16 

29                       1 

7.625 

.3  75 

5C 

17 

29                        1 

7.625 

.375 

51 

20 

32                        1 

11.625 

.3  75 

52 

19 

30                       1 

11.625 

.375 

e-a 

20 

31                        1 

7.625 

.375 

5R 


THE   3R0 


3 

3 
3 
3 

3 
3 

1 
I 
1 
I 
5 
5 
5 
5 
5 
5 
5 
7 
7 
7 
7 
7 
7 
7 
7 
7 

J   1 

1 

1 
1 
1 
9 
9 
7 
7 
9 


2 
2 
2 
2 
6 
6 
6 
6 
6 
6 
6 
3 

a 

9 

a 

9 
3 
3 
3 
3 

2 

2 

2 

2 

10 

10 

3 

3 

10 
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TABLE  XVIII 
DATA  FILE  DD-9  6  3  FOREMAST  continued 


54 

19 

31                        1 

7.625 

.3  75 

55 

18 

24                        1 

7.625 

.3  75 

56 

21 

24                        1 

7.625 

.375 

57 

24 

28                        1 

7.625 

.375 

58 

24 

29                        1 

7.625 

.3  75 

59 

16 

25                        1 

11. 5G 

.5C 

60 

17 

25                        1 

11. 5C 

.  5C 

MEMBERS    61    ThRU    70 

MGCEL    THE    SPQ-9    PLA" 

61 

26 

29                         1 

11.50 

.50 

62 

27 

28                        1 

11.50 

.50 

62 

29 

32                        1 

11.50 

.50 

64 

28 

30                        1 

11.50 

.50 

65 

31 

32                        1 

11.50 

.  5C 

66 

31 

30                        1 

11.50 

.50 

67 

28 

29                        1 

11.50 

.50 

ELEMENTS    63 

*N0    69 

APE    DUMMY    ELEMENTS 

it 

15 

60                         1 

.CC1 

.0001 

69 

15 

60                         1 

.CC1 

.0001 

70 

30 

32                         1 

11. 5C 

.  50 

ENC    THE    3RD 

TIER    CF 

THE    GUAORAPOC    STAR' 

71 

29 

39                         1 

11.75 

.  75 

72 

28 

27                         1 

11.75 

.75 

73 

26 

34                         1 

11.75 

.75 

74 

27 

36                        1 

11.75 

.  75 

75 

26 

39                         1 

7.625 

.375 

76 

27 

37                         1 

7.625 

.3  75 

77 

29 

40                         1 

7.625 

.3  75 

78 

28 

38                         1 

7.625 

.375 

79 

26 

52                         1 

7.625 

.375 

eo 

27 

53                        1 

7.625 

.3  75 

El 

28 

33                        1 

7.625 

.3  75 

62 

29 

33                        1 

7.625 

.375 

83 

■33 

37                        I 

7.625 

.375 

84 

-  ■a 

39                        1 

7.625 

.3  75 

65 

34 

39                         1 

7.625 

.375 

€6 

36 

37                        1 

7.625 

.375 

87 

39 

40                         1 

7.625 

.375 

88 

37 

38                    '    1 

TFCRM 


THE 


9 

10                       1 

9 

10                        1 

9 

10                        1 

9 

10                       1 

9 

10                         1 

9 

10                       1 

9 

10                       1 

ALGNG    WITH 
11 

A    FIXEO    MASS 
12                       1 

11 

12                       1 

11 

12                       1 

11 

12                        1 

11 

12                       1 

11 

12                        1 

11 

12                        1 

0 

0                       1 

0 

0                        1 

11 

12                       1 

4TH    TIER 
1 

2                       1 

1 

2                       1 

1 

2                        1 

1 

2                        1 

13 

14                       1 

13 

14                       1 

13 

14                        1 

13 

14                        1 

13 

14                        1 

13 

14                        1 

13 

14                        1 

13 

14                         1 

13 

14                        1 

13 

14                         1 

13 

14                         1 

13 

14                         1 

13 

14                        1 

13 

14                          1 

85 


TABLE  XIX 


DATA  FILE  DD-9  6  3  FOREMAST  continued 


7.625 

.375 

124 

53 

54                         1 

15 

16                         1 

7.625 

.3  75 

125 

42 

54                        1 

15 

16                         1 

7.625 

.375 

126 

43 

54                         1 

15 

16                         1 

7.625 

.375 

127 

42 

44                         ] 

15 

16                         1 

7.625 

.375 

128 

43 

59                         1 

15 

16                         1 

7.625 

.3  75 

5    ENC    THE    LOMEF 

CJADPAPOD 

START    THE    TOP^ 

AST    MEMBERS 

129 

52 

58                         1 

17 

18                         1 

7.625 

.375 

no 

52 

57                        1 

17 

18                         1 

7.625 

.375 

121 

53 

58                        1 

17 

18                         1 

7.625 

.375 

122 

57 

61                   •      ] 

17 

18                         1 

7.625 

.375 

123 

58 

61                         1 

17 

18                         1 

7.625 

.3  75 

124 

57 

58                       1 

17 

18                         1 

7.625 

.375 

125 

57 

63                         1 

17 

18                         1 

7.625 

.375 

126 

58 

63                         ] 

17 

18                        1 

7.625 

.375 

127 

14 

25                        1 

19 

20                        1 

23. CO 

1.5 

126 

14 

62                        1 

19 

20                        1 

23. CO 

1.5 

125 

62 

61                         ] 

19 

20                        1 

23. CC 

1.5 

140 

61 

63                         1 

L                     19 

20                        1 

22. CO 

1.5 

141 

63 

64                         1 

19 

20                        1 

23. CO 

1.5 

1*2 

64 

65                         1 

21 

22                        1 

5.50 

.50 

1*3 

65 

66                         ] 

I                     21 

22                        1 

5.f0 

.  5C 

144 

66 

67                         ] 

L                     21 

22                        1 

5.50 

.  50 

145 

27 

29                         ] 

L                     13 

14                        1 

7.625 

.3  75 

!    (jPPEP     >ARDAR* 

S    ANC    BRACE 

c 

146 

64 

70                         ] 

L                    23 

24                       1 

14. CC 

.75 

147 

64 

71                        1 

23 

24                         1 

14. CO 

.75 

148 

70 

72                        ] 

L                     25 

26                         1 

1C.CC 

.50 

149 

71 

73                        ] 

L                    25 

26                        1 

1C.CC 

.  50 

150 

72 

68                        ' 

L                    25 

26                        1 

10. CO 

.50 

151 

72 

69                        ) 

L                     25 

26                        1 

10. CO 

.50 

152 

61 

72                        ] 

L                     27 

28                        1 

7.625 

.375 

153 

61 

73 

L                     27 

28                        1 

7.625 

.3  75 

154 

63 

70                    ; 

L                     29 

30                        1 

5.75 
155 

.25 
63 

71                     : 

L                     29 

30                         1 

5.75 
156 

.25 

61 

72 

I                     27 

28                        1 

7.625 

.375 

157 

61 

73 

L                     27 

28                        1 

7.625 

.3  75 

158 

63 

70 

L                     29 

30                        1 

36 


— '.'  "  ~  ~ 


7.625 

.375 

i 

=V-     -_:     _: : 

TIB     CF 

7J-S    CJASa-iPGC 

START 

w 

■5  C 

47 

11.75 

.75 

SC 

37 

45 

11.15 

.75 

SI 

3* 

42 

11.75 

.75 

<2 

36 

43 

11.75 

.  "5 

e-s 

34 

47 

".625 

.3  75 

S4 

36 

45 

7.625 

.3"5 

«e 

4C 

43 

7.625 

.375 

c- 

36 

46 

7.625 

.3  75 

<7 

25 

26 

7.625 

.375 

sa 

25 

27 

7.625 

.375 

c? 

26 

35 

7.625 

•375 

::: 

27 

25 

7.625 

.375 

IC1 

3* 

35 

7.625 

.3  75 

1C2 

35 

36 

7.425 

.375 

::2 

•:  e 

*2 

7.625 

.37  = 

IC4 

35 

43 

7.625 

.375 

1C5 

37 

41 

".625 

.375 

ice 

39 

41 

7.6  25 

.375 

IC1 

*1 

45 

7.625 

.375 

:c3 

-1 

47 

'.625 

.3-5 

< 

"=y=E'£    139 

THRU    123 

"CD5L    "-=     SPG 

-60    s 

l« 

42 

47 

I 

7.625 

.3  "5 

I :: 

-  3 

45 

I 

".625 

.375 

111 

<*Z 

43 

1 

7.625 

.375 

: :; 

45 

47 

1 

7.625 

.375 

113 

*€ 

47 

1 

7.625 

.3  75 

114 

;t 

46 

1 

1  .625 

.3"! 

■  *  5 

tfl 

M 

1 

7.625 

.375 

: :: 

52 

53 

1 

7.625 

.375 

117 

26 

51 

| 

7.625 

.375 

■  :• 

27 

51 

1 

-.625 

.3"5 

■  •  c 

52 

51 

| 

7.425 

.375 

120 

52 

51 

1 

".625 

.375 

121 

44 

52 

1 

7.625 

.375 

122 

=  •3 

59 

1 

7.625 

.3  75 

123 

52 

54- 

1 

>s 


5TH    Tl  l' 

1 

2 

1 

2 

1 

2 

1 

2 

15 

16 

15 

16 

13 

14 

13 

14 

13 

14 

13 

14 

13 

14 

13 

14 

13 

14 

13 

14 

13 

M 

13 

14 

13 

14 

13 

14 

13 

14 

13 

iu 

R"    ALCNG 

15 

ilTfc   i 

16 

15 

16 

i5 

16 

15 

16 

15 

16 

15 

16 

15 

16 

15 

16 

13 

14 

13 

14 

13 

14 

13 

14 

15 

16 

15 

16 

15 

16 

FIXED 


iSS 


TABLE  XXI 
DATA   FILE  DD-96  3  FOREMAST  continued 


5.75 

.25 

159 

63 

71 

1 

29 

30                       1 

5.75 

.25 

i    SPS55    PLATFQPM 

160 

63 

22 

1 

29 

30                       1 

5.75 

.25 

HI 

63 

23 

1 

29 

30                      1 

5.75 

.25 

162 

=     is 

22 

5  C 

23 

1 

29 

30                        1 

-  .  1 3 

.  c  - 

163 

64 

23 

1 

29 

30                       1 

5.75 

.25 

164 

64 

22 

1 

29 

30                         1 

5.75 

.25 

$    LOV.ER    >ARDAR* 

ANO    BRACES 

165 

44 

62 

1 

15 

16                         1 

7.625 

.375 

166 

59 

62 

1 

15 

16                         1 

7.625 

.375 

167 

49 

44 

1 

31 

32                         1 

11.  50 

.50 

168 

59 

50 

1 

31 

32                        1 

11.50 

.50 

16"? 

49 

55 

1 

31 

32                        I 

11.50 

.50 

no 

50 

56 

1 

31 

32                       1 

11. 5C 

.  5C 

171 

14 

49 

1 

33 

34                       1 

5.75 

.25 

172 

14 

50 

1 

33 

34                        1 

5.75 

.25 

<«**** 44 41****4**4*2*44**44**3**** 2****44*  ****  *4*5  **  **4*4  4  *644*  4444  4  4  7 

J6L0CK     I    JOINT    CONSTRAINT 

OATA 

NCJ    CARDS 

FORMAT 

7110 

S       JN                        I> 

IY 

IZ 

I  YY 

IZZ 

J                        1 
1 

1 

1 

1 

I 

\ 

1             I 

2 

1 

1 

1 

I 

1                         1 

* 

1 

1 

1 

1 

1                         1 

4 

1 

1 

1 

1 

1                         1 

5 

0 

0 

0 

0 

0                       0 

6 

c 

0 

0 

0 

0                       0 

7 

0 

0 

0 

0 

0                       0 

8 

0 

0 

0 

0 

0                        0 

9 

c 

0 

0 

0 

0                        0 

10 

0 

0 

0 

0 

0                       0 

11 

0 

0 

0 

0 

0                       0 

12 

c 

0 

0 

0 

0                        0 

13 

0 

0 

0 

0 

0                       0 

14 

0 

0 

0 

0 

0                       0 

15 

1 

1 

1 

1 

1                1 

16 

0 

0 

0 

0 

0                       0 

17 

0 

0 

0 

0 

0                       0 

IS 

0 

0 

0 

0 

0                       0 

19 

0 

0 

0 

0 

0                         0 

20 

0 

0 

0 

0 

0                        0 

21 

c 

0 

0 

0 

0                        0 

22 

0 

0 

0 

0 

0                       0 

23 

0 

0 

0 

0 

0                        0 

24 

0 

0 

0 

0 

0                        0 

25 

0 

0 

0 

0 

0                        0 

26 

0 

0 

0 

0 

0                        0 

27 

0 

0 

0 

0 

0                        0 

28 

0 

0 

0 

0 

0                       0 

29 

0 

0 

0 

0 

0                       0 

20 

c 

0 

0 

0 

0                       0 

31 

0 

0 

0 

0 

0                       0 

22 

0 

0 

0 

0 

0                        0 

23 

0 

0 

0 

0 

0                        0 

34 

0 

0 

0 

0 

0                        0 

i  e 

0 

0 

0 

0 

0                        0 

26 

0 

0 

0 

0 

0                       0 

37 

0 

0 

0 

0 

0                       0 

444 


TABLE  XXII 
DATA  FILE  DD-9  6  3  FOREMAST  continued 


38 

0 

0 

0 

0 

0 

0 

39 

0 

0 

0 

0 

0 

0 

40 

0 

0 

0 

0 

0 

0 

41 

0 

0 

0 

0 

0 

0 

42 

c 

0 

0 

0 

0 

0 

A3 

0 

0 

0 

0 

0 

0 

44 

0 

0 

0 

0 

0 

0 

45 

0 

0 

0 

0 

0 

0 

46 

0 

0 

0 

0 

0 

c 

47 

0 

0 

0 

0 

0 

0 

48 

c 

0 

0 

0 

0 

0 

49 

0 

0 

0 

0 

0 

0 

50 

c 

0 

0 

0 

0 

0 

51 

0 

0 

0 

0 

0 

0 

52 

0 

0 

0 

0 

0 

0 

53 

c 

0 

0 

0 

0 

0 

54 

o 

0 

0 

0 

0 

0 

55 

0 

0 

0 

0 

0 

0 

56 

c 

0 

0 

0 

0 

0 

57 

0 

0 

0 

0 

0 

0 

58 

0 

0 

0 

0 

0 

c 

59 

c 

0 

0 

0 

0 

0 

CO 

I 

1 

1 

1 

1 

1 

s 

0 

0 

0 

0 

0 

0 

c 

0 

0 

0 

0 

0 

63 

0 

0 

0 

0 

0 

0 

64 

0 

0 

0 

0 

0 

0 

65 

0 

0 

0 

0 

0 

0 

66 

0 

0 

0 

0 

0 

0 

C7 

0 

0 

0 

0 

0 

0 

68 

0 

0 

0 

0 

0 

0 

69 

0 

0 

0 

0 

c 

0 

7C 

G 

0 

0 

0 

0 

0 

71 

0 

0 

0 

0 

0 

0 

72 

0 

0 

0 

0 

0 

0 

73 

0 

0 

0 

0 

0 

0 

$4*44*  4441***  4  <*44c*2*********3*********4**  4******5*********6*********  7: 

J3LCCK    J    JOINT 

LCADING 

CATA 

1 

CARD 

FORMAT 

no 

1                (CfIT     IF 

NLC; 

=c   n 

8LCCK 

$       NLJ 
S                       1 
76 

i 

NLJ 

CARDS 

FORMA 

S          JN 

FX 

FY 

FZ 

TX 

TY 

TZ 

«    l-ALYARC    AND 

* : -  = 

AMENf 

5                           1 

1 

1 

1 

1 

1 

! 

49 

-5C. 

-100, 

.0 

0, 

,0 

0. 

0. 

0. 

50 

-50. 

-100 

.0 

0. 

,0 

0. 

G. 

0. 

55 

-50. 

-100 

.0 

0, 

.0 

0. 

0. 

0. 

56 

-50, 

-ICO, 

.0 

0. 

.0 

0. 

0. 

0. 

67 

1 

D. 

-50 

.0 

0. 

.0 

0. 

0. 

c. 

J     SHIP'J    HOTlCh    LOAD 

S          JN 

t                           1 
1 

FX 

1 

FY 

1 

FZ 

1 

TX 

1 

TY 

1 

TZ    1 

0 

.0 

0 

.0 

160, 

1 

.5 

0.0 

0.0 

0.0 

2 

0 

.c 

0, 

.0 

163. 

,  3 

G.O 

0.0 

0.0 

3 

0 

.0 

0 

.0 

96, 

.4 

0.0 

0.0 

0.0 

4 

0, 

>c 

0, 

.0 

137, 

.5 

0.0 

0.0 

0.0 

5 

0 

.0 

0 

.0 

1C8. 

.  2 

0.0 

0.0 

0.0 

6 

0 

.0 

0 

.0 

194, 

.3 

0.0 

0.0 

0.0 

7 

o, 

.0 

o. 

.0 

66. 

.7 

0.0 

0.0 

0.0 

8 

0 

.0 

0 

.0 

164, 

.3 

0.0 

0.0 

G.G 

9 

0 

.c 

0 

.0 

414, 

.9 

0.0 

0.0 

0.0 

10 

0 

.0 

0 

.0 

467. 

.0 

0.0 

G.O 

0.0 

11 

0 

.0 

0 

.0 

381, 

.6 

0.0 

0.0 

0.0 

12 

0. 

.c 

0, 

.0 

381. 

.6 

0.0 

0.0 

0.0 

13 

0 

.0 

0 

.0 

70. 

,  2 

G.O 

0.0 

0.0 

14 

0 

.0 

0 

.0 

325, 

.4 

0.0 

0.0 

0.0 

16 

0, 

.0 

0, 

.0 

565, 

.  7 

0.0 

0.0 

0.0 

17 

0 

.0 

0 

.0 

528, 

.5 

0.0 

0.0 

0.0 

18 

0 

.G 

0 

.0 

4  29 

.1 

0.0 

0.0 

0.0 

19 

0 

.0 

0 

.0 

2  26. 

.9 

0.0 

0.0 

0.0 

39 


TABLE  XXIII 
DATA  FILE  DD-96  3  FOREMAST  continued 


20 

0.0 

0.0 

214.8 

0.0 

0.0 

0.0 

21 

0.0 

0.0 

468.9 

0.0 

0.0 

0.0 

22 

0.0 

0.0 

529.0 

0.0 

0.0 

0.0 

23 

0.0 

0.0 

5  29.0 

0.0 

0.0 

cc 

24 

o.c 

0.0 

91.6 

0.0 

0.0 

0.0 

25 

o.c 

0.0 

4  74.7 

0.0 

0.0 

0.0 

26 

0.0 

0.0 

3234.2 

0.0 

CO 

CO 

27 

0.0 

0.0 

2296.2 

0.0 

0.0 

0.0 

28 

0.0 

0.0 

18  74.0 

0.0 

0.0 

0.0 

29 

0.0 

0.0 

1874.0 

0.0 

0.0 

0.0 

30 

0.0 

0.0 

1459.2 

0.0 

0.0 

0.0 

21 

0.0 

0.0 

1238.5 

0.0 

0.0 

CO 

22 

o.c 

0.0 

14  89.9 

0.0 

0.0 

0.0 

23 

0.0 

0.0 

211.7 

0.0 

0.0 

0.0 

24 

0.0 

0.0 

498.7 

0.0 

0.0 

CO 

25 

o.c 

0.0 

606. 8 

0.0 

0.0 

0.0 

26 

0.0 

0.0 

498.  7 

0.0 

CO 

0.0 

37 

0.0 

0.0 

725.9 

0.0 

0.0 

0.0 

28 

0.0 

0.0 

174.3 

0.0 

CO 

0.0 

29 

0.0 

0.0 

725.9 

0.0 

0.0 

0.0 

40 

o.c 

0.0 

174.3 

0.0 

0.0 

0.0 

41 

o.c 

0.0 

345.5 

0.0 

CO 

0.0 

42 

0.0 

0.0 

1729.0 

0.0 

0.0 

CO 

43 

0.0 

0.0 

1729.0 

0.0 

0.0 

0.0 

44 

0.0 

0.0 

1425.  6 

0.0 

0.0 

0.0 

45 

0.0 

0.0 

1810.4 

0.0 

0.0 

CO 

46 

0.0 

0.0 

1468.1 

0.0 

0.0 

0.0 

47 

0.0 

0.0 

1810.4 

0.0 

0.0 

CO 

46 

0.0 

0.0 

1468.1 

0.0 

0.0 

0.0 

49 

o.c 

0.0 

185.6 

0.0 

0.0 

0.0 

50 

0.0 

0.0 

185.6 

0.0 

CO 

CO 

51 

o.c 

0.0 

540.9 

0.0 

0.0 

0.0 

52 

0.0 

0.0 

965.  3 

0.0 

0.0 

CO 

53 

0.0 

0.0 

783.4 

0.0 

0.0 

0.0 

54 

0.0 

0.0 

464.2 

0.0 

0.0 

0.0 

55 

0.0 

0.0 

113.0 

0.0 

0.0 

CO 

56 

0.0 

0.0 

113.0 

0.0 

0.0 

0.0 

57 

0.0 

0.0 

294.5 

0.0 

0.0 

0.0 

58 

0.0 

0.0 

426.2 

0.0 

CO 

0.0 

59 

0.0 

0.0 

1425.6 

0.0 

0.0 

0.0 

61 

0.0 

0.0 

723.  8 

0.0 

CO 

0.0 

62 

0.0 

0.0 

3  04.6 

0.0 

0.0 

0.0 

63 

0.0 

0.0 

544.  C 

0.0 

0.0 

0.0 

64 

0.0 

0.0 

14C7.7 

0.0 

0.0 

0.0 

65 

0.0 

0.0 

125.0 

0.0 

0.0 

0.0 

66 

0.0 

0.0 

128.8 

0.0 

0.0 

0.0 

67 

0.0 

0.0 

100.3 

0.0 

0.0 

o.c    . 

68 

o.c 

0.0 

60.1 

0.0 

0.0 

0.0 

69 

0.0 

0.0 

60.  1 

0.0 

0.0 

0.0 

70 

0.0 

0.0 

463.2 

0.0 

0.0 

0.0 

71 

o.c 

0.0 

1713.1 

0.0 

CO 

0.0 

72 

0.0 

0.0 

227.0 

0.0 

CO 

CO 

73 

0.0 

0.0 

1081.2 

0.0 

0.0 

CO 

!>«»j<«l«i»4«t*i«44i***as3  1?»*»»«tud«  ****  A**  5*-*  *******  6***  ******  7*** 

S8LCCK    K    FTXEC 

MASS     CATA 

NFMASS    CAPOS 

FORMAT 

IlOtFlCO 

$               (CNI7    IF 

\FMASS=C    IN 

BLOCK 

31 

S          JN                  MA 
<                        | 

ss 

$         PLATFORM    3ELCW    SPQ-5 

18 

32.0 

19 

32.0 

2C 

22.  C 

21 

22.0 

i         SPC-9    AND 

FLATFCPM 

26 

S90.0 

27 

s«o.c 

28 

890  .0 

29 

890.0 

20 

eso.c 

21 

e«=o.c 

22 

890.0 

SPG-60    ANC    PLATFCRP 


90 


42 

«40.0 

43 

<=40.C 

44 

<=40  .0 

45 

<?4O.0 

46 

=  40. C 

47 

S40.C 

4e 

«40.0 

55 

<540  .0 

1         SPS-55   AND 

PLATFCRM 

22 

325.0 

23 

325.0 

64 

225.  C 

f         ANTENNAS 

ANC 

49 

25.0 

fC 

25.0 

55 

12.0 

56 

12.0 

64 

220.0 

67 

20.  C 

68 

26.  C 

69 

26.0 

7C 

120.0 

71 

C72.C 

73 

565. C 

$44*4*4  441*****  44*4?  *44  4*44**-)******  44*4**  4***  444  5  4*4******  6****  *****  7*** 

SBICCK    L    CESIGN    VARIABLE 

INFCRMATIO 

1       C  AP  0  S 

AS     REQ'O            PCRMAT 

3F10.0 

1            (AREA    ANC 

CIMENSICN 

VARIABLES 

GMI  T 

IF     N0VAR1=0» 

S            XMI)    . 

.       .       XA(NDVARl) 

*           13.4 

.tJ 

7.625 

.371 

8.625                 .375 

7.625 

.375 

8.625 

.375 

12.50 

.50 

7.625                  .375 

7.625 

.275 

7.625 

.375 

23.00 

1.50 

9.50                    .50 

14.00 

.75 

10.00 

.50 

7.625 

.375 

7.625                 .375 

11.50 

.50 

5.75 

.25 

$       xAim 

1 

I 

1 

1                           1 

1 

8. CO 

.50 

5.000 

.25  0 

5.000                  .250 

5.00C 

.250 

5. CCO 

.250 

8.00 

.25 

5.000                 .250 

5.0CC 

.250 

5.0CO 

.2  50 

15.00 

.50 

6.00                    .25 

10.00 

.50 

6. CO 

.25 

5.000 

.250 

3.00                    .10 

8.00 

.25 

2. CO 

.  10 

$         XAUdj 
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APPENDIX  3 
PROGRAM  ORGANIZATION 

A.   DESCRIPTION 

The  program  organization  is  layed  out  in  the  following 
flow  charts.   The  main  driver  program  (SADXM)  calls  a 
subdriver  (3AIXSD),  and  the  optimizer  of  the  user's  choice. 
All  changes  required  for  replacement  of  the  optimizer  are 
made  in  SADXM,  thus  allowing  for  easy  testing  of  several 
optimizers  on  the  same  problem. 

SADXSE  may  be  called  from  the  main  fcr  input,  analysis, 
and  cutout.   Printed  output  may  vary  as  the  user  requires. 
A  complete  listing  of  all  subroutines  and  their  functions  is 
ffiven  in  Table  XXIV. 
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TABLE    XXIV 

SUBROUTINE    DIRECTORY 

3T"D"0  DRIVER    PROGRAM~FO"R_USING    THE    ABOVE 

SUBROUTINES.     SADXM    MAY    BE    COUPLED 
TO    OPTIMIZER    OF    USER'S    CHOICE. 

SADXSD  SU5DRIVER    PROGRAM    FOR    COUPLING 

SACX    ROUTINES    TO    SADXM 

SADX01  THIS    ROUTINE    READS    AND    PRINTS    INPUT    DATA 

AND    ORGANIZES    PSEUDO-DYNAMIC    STORAGE 
ALLOCATION 

SADX02  3UILDS    VECTORS    JC    AND    UK    FOR    FINITE    ELEMENT 

STRUCTURAL    ANALYSIS 

SADX03  BUILDS    THE    12x12    ELEMENT    STIFFNESS    MATRIX 

SADX05  SUPERIMPOSES    THE    ELEMENT    STIFFNESS  MATRIX    BK 

(OR    ELEMENT     MASS    MATRIX    EM)     ON    THE  COMPACTED 

GLOBAL    STIFFNESS    MATRIX    AK     (OR    THE  GLOBAL 
COMPACTED    MASS    MATRIX    AM) 

SADX06  BUILDS    GLOBAL    LUMPED    MASS    MATRIX 

5ADX07  LU    DECOMPOSES    SYMMETRIC, POSITIVE-DEFINITE, 

SPARCE  MATRICESrTHE  UPPER  TRIANGLE  0?  WHICH 
IS  STORED  IN  MATRIX  AK  (OR  AM)  WITH  LEADING 
ZEROES    NOT    STORED 

SADX08  FORWARD    AND    BACK    SUBSTITUTES    TO    YIELD    A 

SOLUTION    A    SET    Or    LINEAR    EQUATIONS 
(DECOMPOSED     SY    SADX07    DR    EQUIVALENT) 

SADX09  PRINTS     ALL    JOINT    DISPLACEMENTS    FOR    EACH 

LOAD    CONDITION    OF    A    FINITE    ELSM.     STRUCTURE 

SADX11  ROUTINE   TO    ORGANIZE    ANALYSIS 

SADX15  CALCULATES    VALUES    FOR    ALL    DESIGN    AND 

3EHAVI0RIAL     CONSTRAINTS    AS    DEFINED    BY 
"SADX"    PR03RAM 

SADX16  CALCULATES    STRESS    IN    TRUSS    ELEMENT    LNO 

UNDER    LOAD    CONDITION    JJ 

SADX17  PRINTS     STRESSES    AND/OR    FORCES    FOR 

TRUSS    ELEMENTS 

SADX19  ADDS    ELEMENT    MASS    MATS  IX    A  A    OF    ELEMENT 

LNO    TO    GLOBAL    MASS    MATRIX    AM    TO    BUILD 
THE    LUMPED    MASS    MATRIX 


SADX23 


CALCULATES  WEIGHT  OF  TRUSS/FRAME  STRUCTURE 
OR  CALCULATE  WEIGHT  DF  INDIVIDUAL  MEMBERS 


SADX36       CALCULATES  (  XEIG-T*AM*XEIG)  FOR 

GRADIENT  CALCULATIONS  IN  FREQUENCY 
CONSTRAINTS 

SADX37       CALCULATES  EIGENVALUE  GRADIENT  INFORMATION 
IN  FINITE  ELEMENT  STRUCTURAL  ANALYSIS  AND 
DESIGN 
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SADX46       READS  INPUT  INFORMATION  FOR  TPUSS 
ELEMENTS 

SADX47       TRANSFORMS  THE  ELEMENT  STIFFNESS  MATRIX 

SK  (OR  ELEMENT  MASS  MATRIX  EM)  FROM  LOCAL 
TO  GLOBAL  COORDINATES 

SADX49       SOLVES  REAL  EIGENVALUE  PROBLEMS  USING 
THE  SUESPACE  ITERATION  METHOD 

SADX50       BUILDS  INITIAL  SET  OF  BASIS  VECTORS 
FOR  EIGENSOLUTION  BY  REDUCED  BASIS 
METHOD 

3ADX53       PRINTS  MSM3ER  INFORMATION  FOR  TRUSS 
ELEMENTS 

SADX62       PRINTS  NEIG  EIGENVALUES  STORED  IN 
EIGVAL,  AND  THEIR  CORRESPONDING 
EIGENVECTORS  STORED  IN  XEIG 

SADX71       PRINTS  G  VECTOR  OF  CONSTRAINTS 

SADX72       READS  IN  FRAME  ELEMENT  INPUT  DATA 

SADX78       BUILDS  3x3  TRANSFORMATION  ARRAY  TRFORM 
FOR  TRANSFORMING  FROM  LOCAL  TO  GLOBAL 
COORDINATES 

SADX80       CALLS  SADX03  TO  3UILD  THE  ELEMENT  STIFFNESS 

MATRIX;  CALLS  SADX47  TO  TRANSFORM  THE  MATRIX; 
AND  CALLS  3ADX05  TO  3UILD  THE  REDUCED  GLOBAL 
STIFFNESS  MATRIX 


SADX81       CALCULATES  STRESS  IN  FRAME  ELEMENT  LNO  UNDER 
LOAD  CONDITION  JJ  (IF  JJ=0  STRESSES  CALCUL- 
ATED FOR  ALL  LOAD  CONDITIONS) 


SADX32       READS  INPUT  DATA  FOR  FRAME  ELEMENTS  WITH 
SPECIFIED  SECTION  TYPES 

SADX8  3       PRINTS  STRESSES  AND/OR  FORCES  FOR  FRAME 
ELEMENTS 

SADX84       PRINTS  MEMBER  INFORMATION  FOR  FRAME  ELEMENTS 

SADX85       CALCULATES  SECTION  PROPERTIES  FOR  FRAME 

ELEMENTS  OF  A  SECTION  TYPE  GIVEN  BY  LSECT 

SADX86       CALLS  EITHER  SADX06  TO  5UILD  THE  LUMPED  MASS 
MATRIX  OR  BUILDS  THE  CONSISTENT  MASS  MATRIX 
BY  CALLING  SADX87  TO  BUILD  THE  ELEMENT  MASS 

SS 
ED 
GLOBAL  "MASS  MATRIX 

SADX87       CALLS  BUILDS  THE  ELEMENT  CONSISTENT  MASS 
MATRIX 
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SADX88 

SADX89 
SADX90 

SADX91 

SETIME 

GSTIME 


CONVEETS    UNFORMATTED    DATA    TO    FORMATTED    DATA 
IN    FIELDS    OF     10,     EACH    RIGHT    JUSTIFIED    AND 
ACCEPTS    ^COMMENT    CARDS    IN    DATA 

SOLVES    EIGENVALUE    PROBLEM    |A-ALAMBD*Bj     7=0 

SOLVES    EIGENVALUE    PROBLEM 

SOLVES    EIGENVALUE    PROBLEM 

STARTS    EXECUTION    TIMER 
(NONiaSL     LIBRARY) 

STOPS    EXECUTION    TIMER 
(NONIMSL    LIBRARY) 
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TABLE  XXV 
PROGRAM  BLOCK  DIAGRAM 
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